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I.  INTRODUCTION 

The  NOVA-2  (Nuclear  Overpressure  Vulnerability  Analysis,  Version  2) 
computer  code,  documented  in  Reference  1,  was  developed  for  use  in  analyzing 
the  dynamic  response  of  aircraft  structure  to  transient  pressure  loads 
associated  with  the  blast  wave  from  a  nuclear  explosion.  Since  then  a  number 
of  changes  have  been  made  to  the  code  to  increase  its  versatility.  The 
purpose  of  this  report  is  to  compile  these  changes  so  that  they  will  be  more 
readily  available  to  the  user. 

The  first  major  change  to  the  NOVA-2  code  (N0VA-2L,  Refs.  2  and  3)  was 
made  to  provide  the  analyst  with  a  versatile  set  of  pressure  loading  routines 
so  analytical  predictions  could  be  compared  with  experimental  tests  such  as 
those  conducted  in  the  Sandia  Laboratories’  Thunderpipe  shock  tube  facility 
(Ref.  4). 

The  pressure  loading  routines  in  the  N0VA-2L  code  are  analytical  in  form 
and  require  a  certain  amount  of  manual  processing  of  the  test  data.  In  order 
to  eliminate  most  of  this  time-consuming  work,  changes  were  made  to  allow  the 
analyst  to  input  digitized  experimental  pressure  data  directly  from  magnetic 
tape,  and  the  N0VA-2LT  code  (Ref.  5)  was  created.  All  of  the  pressure  load 
options  from  NOVA-21  are  contained  in  N0VA-2LT  unchanged.  The  changes 
necessary  to  transform  NOVA-2  to  N0VA-2LT  are  listed  in  Appendix  A. 

Similar  changes  are  documented  in  Reference  6  for  the  N0VA-2S  code, 
a  stiffened  panel  extension  to  the  basic  NOVA-2  code.  This  special  version 
is  called  N0VA-2LTS  where  LT  refers  to  the  fact  that  all  of  the  nuclear  blast 
and  aeronautical  loading  routines  have  been  replaced  by  user-generated  functions 
or  digitized  test  data  on  tape. 

Additional  changes  were  made  to  the  NOVA-2  (N0VA-2LT)  and  N0VA-2S 
(N0VA-2LTS)  versions  to  include  a  pinned/sl iding-pinned  end  constraint. 

N0VA-2LTS  was  also  revised  to  include  free-free  and  clamped-free  edge  con¬ 
ditions,  and  elastic  springs  along  the  edges  (Ref.  7).  These  changes  are 
documented  in  Appendixes  B  and  C. 

Section  II  describes  the  pressure  load  options  available  to  specify  the 
loading  applied  to  the  structural  element  model. 
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Section  III  covers  program  description  and  operation,  including  rib 
element  end  constraints,  stiffened  panel  edge  fixity,  and  input  data 
instructions. 

Section  IV  describes  a  method  for  approximating  the  curvature  of  a 
noncylindrical  panel  using  geometric  radial  imperfections. 
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II.  DYNAMIC  LOAD  OPTIONS 

Four  transient  pressure  models  and  one  static  pressure  model  are  incor¬ 
porated  into  the  N0VA-2LT  and  N0VA-2LTS  codes.  Two  of  the  transient  models 
are  uniform  spatially,  i.e.,  the  pressure  functions  vary  with  time  only.  The 
static  pressure  model,  of  course,  is  independent  of  time  and  is  represented 
by  specifying  either  a  positive  or  negative  pressure,  P  ,  to  represent  the 
direction  and  magnitude  of  the  loading.  This  is  consistent  with  NOVA-2,  where 
a  positive  pressure  is  taken  to  be  acting  inward.  A  combination  static-dynamic 
response  calculation  can  be  made  where  the  static  solution  is  determined  first, 
and  then  followed  by  the  dynamic  load  and  resulting  response. 

For  cases  where  the  code  is  being  used  in  the  iterative  modes,  range  is 
replaced  by  a  range  factor  parameter  which  is  normalized  to  1.0  initially. 

For  subsequent  trials  this  parameter  acts  like  a  range;  but  it  is  actually 
the  inverse  of  the  factor  which  adjusts  the  pressure  loading.  Only  the 
magnitude  of  the  pressures  is  affected,  not  the  characteristic  times. 

1.  DYNAMIC  LOAD  OPTION  1— SPATIALLY  UNIFORM  LOAD  WITH  AN  ANALYTICAL  DECAY 
FUNCTION. 

The  first  dynamic  load  option  consists  of  a  uniform  pressure  (Px)  with  an 
analytical  decay  function  designed  to  approximate  the  diffraction  and  drag 
phases  of  a  blast  loading.  This  function  (see  Figure  1)  represents  a  combin¬ 
ation  linear  and  exponential  decay: 

Pi  (t)  =  Pi  (l  '  £-)  (t  <  t') 
i 

/  t_  \n  -at/t 

PII(t)  =  Pol1  '  to  )  e  °  (f  it  <  t0) 

Pm(t)  =  0  (t<tQ)  (1) 

The  second  function  of  Equation  1  is  used  for  time  greater  than  or  equal  to  t' . 
By  specifying  t1  -  0,  the  special  loading  cases  indicated  in  Table  1  can  easily 
be  generated,  where  I  is  the  impulse  and  At  is  the  integration  time  interval. 
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Figure  1.  Analytical  pressure  time  history 
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TABLE  1. 

LOADING  PARAMETERS 

FOR  SPECIAL 

CASES 

Step 

Triangular 

Exponential 

Parameter 

Function 

Decay 

Impul se 

Decay 

Pi 

0 

0 

0 

0 

P 

P0 

P 

21 

P 

At 

10 

.20 

t 

0 

1.0  x  10 

t 

0 

1.0  x  10 

lo 

t1 

0 

0 

0 

0 

a 

0 

0 

0 

a 

n 

0 

1 

1 

0 

2.  DYNAMIC  LOAD  OPTION  2--SPATIALLY  UNIFORM  DISCRETE  LOADS. 

Like  the  previous  option,  loading  option  2  is  appropriate  for  uniformly 
applied  loads  without  consideration  of  engulfment.  With  this  option,  discrete 
values  of  pressures  can  be  specified  at  a  set  of  times  beginning  at  zero. 

For  other  times,  linear  interpolation  is  used  except  after  the  last  time  in 
the  table,  where  a  pressure  equal  to  the  last  value  is  assumed. 

3.  DYNAMIC  LOAD  OPTION  3--SPATIALLY  NONUNIFORM  DISCRETE  LOADS. 

The  third  load  option  permits  a  spatially  nonuniform,  discrete  load 
description.  As  in  option  2  the  transient  load  is  specified  at  user-selected 
times. 

For  beam  element  analysis,  the  spatial  variation  is  accomplished  by 
initially  selecting  pressure  measurement  locations  along  the  beam  element. 

These  positions  are  indicated  in  terms  of  mass  point  positions,  i.e.,  position 
3.67  would  indicate  a  measurement  position  2/3  of  the  distance  between  mass 
numbers  3  and  4  in  the  beam  element  model.  Each  measurement  station  then  has 
its  own  time  history. 

A  rectangular  mesh  of  measurement  positions  must  be  used  for  two- 
dimensional,  panel  elements.  This  mesh  need  not  be  the  same  as  the  rect¬ 
angular  integration  grid  used  in  determining  the  structural  response.  Measure¬ 
ment  positions  are  specified  in  terms  of  x,  y  coordinates.  Again,  each  mesh 
point  station  has  its  own  unique  time  history. 
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Interpolations  for  mass  points  or  integration  points  located  between 
measurement  stations  are  made  linearly.  Linear  extrapolation  is  used  for  points 
outside  the  mesh. 

Linear  interpolation  is  also  used  between  break  points  in  the  temporal 
description  except  for  times  greater  than  the  last  table  point,  where  a 
pressure  equal  to  the  last  value  is  assumed.  At  least  one  point  in  the  field 
must  be  intercepted  at  time  zero. 

4.  DYNAMIC  LOAD  OPTION  4--EXPERIMENTAL  TEST  PRESSURE--TIME  HISTORY  FROM 
DIGITIZED  DATA  TAPE. 

The  fourth  load  option  provides  for  user  supplied  pressure-time  histories 
on  magnetic  tape.  After  the  experimental  test  data  have  been  filtered  and 
digitized,  it  is  stored  on  data  tapes  in  a  specific  format  which  is  discussed 
in  Section  III.  Each  tape  contains  multiple  files  where  each  file  corresponds 
to  a  particular  pressure  gage.  This,  of  course,  precludes  the  possibility 
of  iterating  for  critical  response,  because  only  one  set  of  pressure  data  is 
available.  The  data  tape  is  mounted  prior  to  program  execution,  and  the  data 
are  transferred  to  large  core  memory  (available  on  computers  such  as  the 
Control  Data  Corporation  7600  and  CYBER  176)  before  the  structural  response 
analysis  begins. 

During  execution  both  temporal  and  spatial  variation  in  pressure  are 
accomplished  through  linear  interpolation,  though  the  spatial  variation  is 
limited  to  only  one  coordinate  direction.  This  means  that  for  panels  the 
variation  can  either  be  circumferential  or  axial,  but  not  both.  In  case  of  a 
uniform  spatial  distribution,  only  one  pressure  time  history  is  specified  and 
spatial  interpolation  is  not  attempted.  The  temporal  variation  is  limited 
to  the  extent  that  the  times  must  be  regular  intervals  apart  and  the  total 
number  of  samplings  be  within  the  dimensional  limits  in  the  program. 

The  measurement  positions  are  specified  in  the  same  manner  as  described 
in  Section  II  3.  They  must  be  coordinated  with  the  tape  data  on  a  gage-to-gage 
basis  so  that  the  right  data  get  applied  to  the  appropriate  place  on  the  beam 
or  panel.  One  set  of  input  data  specifies  which  tape  file  corresponds  to  which 
measurement  station.  It  is  also  possible  to  exclude  an  entire  set  of  data  if 
that  gage  is  later  found  to  have  been  defective.  The  measurement  stations  are 
assumed  to  be  ordered  consecutively  from  the  beginning  of  the  beam  model,  or 
from  the  coordinate  origin  for  panel  models. 
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In  general,  the  time  history  contained  on  tape  will  include  times  prior 
to  the  shock  arrival  at  the  gage  nearest  the  blast.  Therefore,  it  is 
necessary  to  scan  a  plot  of  the  pressure-time  history  data  obtained  from  that 
gage  to  determine  the  start  time  for  input  of  pressures  into  the  N0VA-2LT(S) 
code. 

The  program  also  provides  for  skipping  over  portions  of  the  pressure-time 
history.  The  input  parameter  (SKIP)  indicates  the  frequency  at  which  the  tape 
data  is  sampled.  For  example,  SKIP=3.0  would  mean  that  every  third  time  point 
following  shock  arrival  (at  the  nearest  gage)  will  be  used  to  describe  the  time 
history.  Using  less  data  should  mean  a  somewhat  faster  running  program,  but 
this  also  sacrifices  some  accuracy.  It  is  the  user's  responsibility  to 
determine  how  much  accuracy  is  required  for  the  pressure-time  history.  It  is 
also  important  that  the  program  integration  stop  time  (TSTOP)  selected  does  not 
exceed  the  data  available  on  the  tape. 
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I  III.  PROGRAM  DESCRIPTION  AND  OPERATION 

The  N0VA-2LT  and  N0VA-2LTS  codes  represent  simplified  versions  of  NOVA-2 
|  since  52  subroutines  have  been  eliminated.  However,  loader  segmentation  is 

still  used  to  optimize  computer  core  requirements  and  flexibility.  If  seg¬ 
mentation  is  not  used,  the  user  must  select  the  necessary  routines  from 
Table  2,  depending  on  which  subprogram  is  being  exercised.  The  user  should 
i  |  realize  that  large  core  memory  (LCM)  is  required  when  using  experimental  test 

:  data  tapes  (load  option  4),  thus  restricting  code  use  to  computers  such  as 

the  Control  Data  Corp  7600  and  CYBER  176  systems. 

The  amount  of  central  memory  required  can  be  reduced  by  separating  the 
code  (N0VA-2LT  or  N0VA-2LTS)  into  two  segmented  codes:  a  beam  element  code 
(DEPROB)  and  a  panel  element  code  (DEPROP). 

1.  BEAM  AND  PANEL  SUBPROGRAMS. 

Since  the  NOVA-2  and  N0VA-2S  program  libraries  are  maintained  at  the  Air 
j  Force  Weapons  Laboratory  (AFWL)  using  the  Control  Data  Corp  6000/7000  UPDATE 

j  program,  it  is  relatively  simple  to  assemble  only  the  subroutines  needed  to 

J  create  the  desired  code  (DEPROB  or  DEPROP)  along  with  the  changes  necessary 

to  transform  the  basic  codes  to  the  specialized  LT  versions  discussed  in 
Section  II.  These  subroutines  are  categorized  in  Table  2.  Note  that  certain 
subroutines  are  common  to  both  codes. 

2.  RIB  BUCKLING  CHANGES  FOR  PINNED  BOUNDARY  CONDITIONS. 

The  N0VA-2LT  (N0VA-2LTS)  codes  have  been  modified  to  include  pinned/ 
si iding-pinned  end  constraints  for  the  rib  option  (KTYPE=10).  Provision  is 
also  made  for  including  the  mass  of  the  piston  and  loading  block  in  the  ana¬ 
lytical  model  to  simulate  the  test  columns  in  the  Boeing  Phase  II  STRESN0 
program  (Ref.  4). 

The  rib  buckling  changes,  documented  in  Appendix  B,  include  a  statement. 
"SM(1)=XXX/DELTS(1)M ,  where  "XXX"  represents  the  mass  in  lb-s2/in.  Since 
this  first  mass  point  actually  represents  a  piston  and  loading  block  in  the 
model,  the  link  between  the  first  end  constraint  at  boundary  coordinates 
(VI,  Wl)  and  the  first  mass  point  at  coordinates  [V(l),  W(l)]  is  a  pseudo  link. 
Initially,  this  link  must  be  horizontal  as  shown  in  Figure  2. 
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TABLE  2.  SUBPROGRAMS  FOR  NOVA-2LT  (NOVA-2LTS) 


COMMON 

SUBROUTINES 

DEPROP 

DEPROB 

NOVA 

DEPROP 

DEPROB 

SEC 

BOLT 

COMP1 

RITER 

aDERVl 

C0MP2 

CSETUP 

DERV2 

COMSET 

PINIT 

DSET1 

CYCLE 

SOLVE 

DSET2 

DAB 

INTL 

DSET3 

DEFORM 

PRESS 

DTSTEP 

DPUR 

HIM 

EQUILP 

LEGEND 

EQUILX 

LIST1 

FB 

LIST2 

FBCTL 

aMATXIN 

FBSET 

RELAXP 

FINAL 

SIGMA 

FSOL 

aSIGMAB 

PRINT1 

aSTIFF 

READ1 

RESD 

RESET 

RLAXB 

RLAXF 

SLAY 

STRESS 

STRESX 

STRN1 

STRN2 

STSET 

TSTEP 

VCS 


aRequired  for  stiffened  panel  version  only. 
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CLAMPED  /  SLIDING  -CLAMPED 


PINNED  SLIDING  -  PINNED 


Figure  2.  Rib  buckling  model. 


1 


AFWL-TR-79-160 

Initial  geometric  imperfections  (in  the  w  direction)  are  required  in 
the  rib  to  induce  buckling.  If  the  default  option  is  used,  a  1-cosine  mode 
shape  is  applied  to  the  clamped  case,  and  a  sine  mode  shape  to  the  pinned  case. 
When  exercising  this  option,  the  input  parameter  AMP  (Group  21,  p.  277,  Ref.  1) 
represents  the  magnitude  of  the  initial  imperfection  at  the  center  of  the  rib. 

The  other  method  of  specifying  an  imperfection  utilizes  the  V(I),  W(I)  mass 
point  coordinate  of  the  rib  (Group  6,  p.  274  of  Ref.  1). 

3.  STIFFENED  PANEL  CHANGES  FOR  FREE  BOUNDARY  CONDITIONS  AND  DISCRETE 
LINEAR  SPRINGS. 

The  N0VA-2S  (N0VA-2LTS)  codes  have  been  modified  to  allow  panels  to  be 
modeled  with  free  edge  boundaries  and  discrete  linear  springs  at  arbitrary 
positions.  These  options  were  added  to  allow  the  user  to  analyze  structures 
such  as  the  center  bomb  bay  doors  of  the  B-52  aircraft.  The  center  door  is 
nearly  clamped  along  one  edge  and  free  along  the  other  three  edges  except  at 
the  two  corners  which  connect  to  the  adjacent  doors.  These  corner  connections 
can  be  simulated  by  springs  representing  the  compliance  of  the  adjacent  doors. 

The  springs  must  be  located  (in  any  of  the  three  orthogonal  coordinate 
directions)  at  the  spatial  integration  points  Yj.3k  defined  for  the  unstiffened 
panel  either  along  the  boundary  edges  or  within  the  panel's  interior.  It  is 
assumed  that  the  coordinate  direction  of  the  spring  forces  are  deformation 
dependent,  that  is,  the  spring  force  always  remains  directed  radially  (w), 
tangentially  (v),  or  axially  (u)  as  the  panel  deforms.  Thus,  the  following 
expression  is  added  to  the  equaton  of  motion  in  Equation  13  of  Reference  6: 


(2) 


where 

K*  =  spring  stiffness  in  w-direction  (lbs/in) 

Sw  =  total  number  of  springs  in  w-direction 

The  other  nomenclature  is  defined  in  References  1  and  6.  Similar  additions 
have  been  made  for  the  u  and  v  equations  of  motion  using  Ku  and  Kv  over  Su  and 
Sv,  respectively.  In  the  future,  rotational  springs  should  be  included  in  a 
similar  manner,  so  that  boundary  conditions  between  clamped  and  simply 
supported  can  be  simulated. 
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The  boundary  combinations  for  the  y  and  p  directions  for  the  w-displacement 
have  been  extended  to  include  the  free-free  and  clamped-free  conditions.  The 
w-displacement  functions  for  the  free-free  condition  in  the  y  and  (3  directions 
are  given  by 


J»w 

m 


Ay  Ay 

=  cosh  -  +  cos  -  -  a  sinh 

7i  7i  m 


Ay  Ay 

m*  m1 

—  +  sin  — 
n  n 


„w  -  i 

=  1  '  l'  ^n  =  2 


<-2 


m  =  3,  4,  5,  .  .  . 


(3) 


Yn0 

=  cosh  +  cos - a  l  sinh 

7T  71  n 


where 


V  +  .  V 
__+  S1n_ 

n  =  3,  4,  5  .  .  . 


A  =  the  roots  of  cosA  coshA  =  1 

_  coshA  -  cosA 
“  sinhA  -  sinA 

The  w-displacement  functions  for  the  clamped-free  condition  in  the  A  and  p 
directions  are  given  by 


Ay  Ay 

„w  _  .  i  •  u 

4»-  =  cosh  -  -  cos - a  \sinh 

m  7t  7X  m 


Ay  Ay 

-  sin  J21 

71  71 

m  =  1,  2,  3, 


(4) 


AS  AS 

<bw  *  cosh - cos  — —  -  a  \sinh 

n  7i  7t  n 


where 


V.sin  V 

7i  n 

n  =  1,  2,  3, 


A  =  the  roots  of  cosA  coshA  =  -1 


a  - 


_  coshA  +  cosA 
sinhA  +  sinA 


For  flat  panels,  these  changes  produce  a  good  approximation  of  the  actual 
boundary  condition,  but  the  bomb  bay  doors  are  curved  and  additional  boundary 
conditions  in  the  u  and  v  directions  are  required.  Presently,  the  u  and  v 
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boundary  conditions  represent  a  held  inplane  condition  along  all  edges.  Held- 
free  and  free-free  boundary  conditions  are  needed  and  should  be  included  in  future 
modificatons  to  the  code. 

The  changes  required  to  extend  N0VA-2S  (N0VA-2SLTS)  to  include  free-free 
and  clamped- free  edge  conditions  and  elastic  springs  are  presented  in 
Appendix  C. 

4.  PROGRAM  SEGMENTATION  AND  OPERATION. 

The  recommended  segmentation  tree  structure  is  presented  in  Figure  3  for 
DEPROP,  and  in  Figure  4  for  DEPROB.  The  common  blocks,  which  must  be  designated 
GLOBAL  and  SAVE  in  the  segmentation  loader  directives,  are  listed  in  Table  3 
along  with  the  routines  to  which  they  are  assigned.  The  subroutines  SIGMA 
and  SIGMAB,  and  the  associated  common  block  CBLANK  are  required  only  for  the 
elastic-plastic  response  option  of  DEPROP.  This  option  corresponds  to  selecting 
NDERV=2  in  Group  5,  page  115  of  Reference  6.  Excluding  them  for  the  elastic 
response  option  (NDERV=1)  decreases  the  amount  of  memory  required  by  23,000g 
cells. 

The  loader  segmentation  directives  for  DEPROP  and  DEPROB  are  listed  in 
Tables  4  and  5  respectively.  Note  that  the  directives  in  Table  4  are  applicable 
to  an  elastic-plastic  panel  solution  since  they  contain  SIGMA,  SIGMAB,  and 
CBLANK. 

The  Fortran  (FTN)  compiler  has  been  used  with  the  NOS/BE  Version  1.2 
Operating  System  to  compile  the  codes  on  the  CYBER  176  computer.  Common  blocks 
are  given  in  Table  6.  In  order  to  reduce  the  LCM  required  and  expedite  running 
time,  the  user  should  always  set  the  P  and  Q  arrays  to  the  amount  of  storage 
needed  for  the  data  tape,  Load  Option  4.  During  operation,  pressure  data  that 
cannot  be  stored  in  P  are  automatically  stored  in  Q.  Check  the  MAXD  and  MAXD2 
dimensions  in  subroutines  PINIT  and  PRESS.  Be  sure  there  is  enough  storage 
capability  for  the  number  of  pressures  selected,  based  on  the  shock  arrival 
time  to  the  first  gage,  the  TSTOP  used,  and  the  time  interval  resulting  from 
the  SKIP  selected. 

The  program  can  currently  handle  11914  pressures  for  each  of  22  channels. 

For  a  sampling  interval  of  10  ps,  119  ms  of  pressure-time  history  can  be 
stored  in  the  program. 
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LEVEL  2 


LEVEL  1 


INTI 

SOLVE 

PRESS 

CSETUP 

YCS 

STRN1 

SLAY 

C0MP2 

COMF1 


STSET 

STRESX 

RLAXB 

RESET 

RESD 

FSOL 

FBSET 

EQUILX 

OPUR 

DAB 


FINAL 

STRESS 

TSTEP 

EQUILP 

READ! 

CYCLE 

CQMSET 
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TA8LE  3.  COMMON  BLOCKS  GLOBAL  AND  SAVE  IN  SEGMENTATION 


Common  block  ownership 
for  segmentation 


Routine 

Common 

Blocks 

'nova 

CNOVA 

CLOAD 

CBLK1 

C0M1 

COM2 

DEPROB 

BLK2 

BLK3 

BLK6 

DEPROP 

CBLK2 

CBLK10 

CBLK3 

CBLK11 

CBLK4 

CBLK13 

CBLK5 

bCBLK15 

CBLK7 

bCBLK17 

bCBLK8 

CBLANK 

CBLK9 

DERV1 

CBLK6 

bCBLK16 

RELAXP 

CBLK12 

aThe  SAVE  designation  is  not  necessary  for  common  blocks  in  the  root  segment. 
Required  for  stiffened  panel  version  (N0VA-2S)  only. 
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TABLE  4.  DEPROP  SEGMENTATION  LOADER  DIRECTIVES 

TREE 

NOVA 

NOVA 

INCLUDE 

PINIT,  SEC 

LEVEL 

TREE 

RITER 

TREE 

DEPROP 

LEVEL 

TREE 

DSETI 

DSETI 

INCLUDE 

DSET2,  LEGEND,  DISTEP 

TREE 

DSET3 

DSET3 

INCLUOE 

BOLT 

aTREE 

aSTIFF 

aSTIFF 

INCLUOE 

aMATXIN 

TREE 

RELAXP 

TREE 

aDERVl 

aDERVl 

INCLUDE 

HIM,  DERV2 ,  LIST1,  LIST2 ,  SIGMA,  aSIGMAB 

LEVEL 

TREE 

CSETUP 

CSETUP 

INCLUDE 

INTI 

TREE 

PRESS 

TREE 

SOLVE 

NOVA 

GLOBAL 

CNOVA,  CLOAD,  CBLK1,  COM1,  COM2 

DEPROP 

GLOBAL 

CBLK2,  CBLK3 ,  CBLK4,  CBLK5 ,  CBLK7 ,  aCBLK8, 

CBLK9 , 

CBLK10 ,  CBLK11, 

CBLK13 ,  aCBLKl 5 ,  aCBLK17 ,  CBLANK-SAVE 

RELAX? 

GLOBAL 

CBLK12-SAVE 

aDERVl 

GLOBAL 

CBLK6,  aC8LK16-SAVE 

END 

aRequired  for  stiffened  panel  version  (N0VA-2S)  only;  otherwise  replace 
DERVI  with  HIM  and  relocate  CBLK6  to  DEPROP  GLOBAL  declaration. 
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TABLE  5.  DEPROB  SEGMENTATION  LOADER  DIRECTIVES 


TREE 

NOVA 

NOVA 

INCLUDE 

PINIT,  SEC 

LEVEL 

TREE 

RITER 

TREE 

DEPROB 

DEPROB 

INCLUDE 

DEFORM,  FB,  FBCTL,  PRINT1,  RLAXF ,  STRN2 

LEVEL 

TREE 

COMP1 

COMF1 

INCLUDE 

C0MP2,  SLAY,  STRN1,  VCS 

-  -TREE  *  ~  " 

DAB  ... 

DAB 

INCLUDE 

DPUR ,  EQUILX ,  FBSET,  FSOL,  RESD,  RESET, 

RELAXB,  STRESX ,  STSET 

TREE 

CYCLE 

CYCLE 

INCLUDE 

EQUIP,  STRESS,  FINAL 

TREE 

COMSET 

COMSET 

INCLUDE 

READ1,  TSTEP 

LEVEL 

TREE 

CSETUP 

CSETUP 

INCLUDE 

INTI 

TREE 

SOLVE 

TREE 

PRESS 

NOVA 

GLOBAL 

CNOVA,  CLOAD,  CBLK1,  COM1,  COM2 

DEPROB 

GLOBAL 

BLK2,  BLK3,  BLK6-SAVE 

END 

NOVA 
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TABLE  6. 

LEVEL  2  VARIABLES 

FOR  LARGE  CORE 

MEMORY  (LCM) 

COMMON 

8L0CKS 

ROUTINES 

VARIABLES 

DIMENSION- 

LENGTH 

FTN 

LIMIT 

LCM 

LIMIT 

C0M1 

PINIT, 

P(NGAGE ,MAXD2) 

22  x  5957 

131,071 

PRESS 

=  131,054 

COM2 

PINIT, 

Q(NGAGE , 

22  x  5957 

131,071 

360,000 

PRESS 

MAXD-MAXD2) 

=  131,054 

CBLANK 

SIGMA 

ALTT(1156),  .  . 

1156  x  14 

131,071 

TTNU( 1156) 

=  16,184 

NOTES:  1.  All  dimensions  in  table  are  in  decimal. 


2.  Common  block  variables  in  CBLANK  can  be  located  in  small  core 
memory  (SCM)  on  the  Control  Data  Corp  6600  computer  with  a  net 
increase  in  core  of  23,000  cells.  Subroutine  SIGMA  is  not 
required  for  beams  or  elastic  panel  response. 

3.  Common  block  variable  P  and  Q  are  required  only  for  load  option 

4;  otherwise,  they  can  be  dimensioned  at  (1,1)  in  either  SCM  or  LCM. 

4.  NGAGE  represents  the  number  of  experimental  channels  of  pressure 
data  on  the  tape.  MAX0  is  the  total  number  of  pressures  to  be 
stored  in  memory.  MAXD2  is  the  number  to  be  stored  in  P;  the 
rest  spill  over  into  Q.  MAXD  can  be  calculated  as  follows: 

MAXD  =  2  +  (TST0P-T1M1)/(SKIP  •  tape  time  interval). 

5.  For  load  option  4,  if  MAXD2  or  MAXD  is  changed,  then  the  two 
statements  beginning  at  statement  number  6000  in  PINIT  need 
to  be  changed. 

6.  Current  limits  on  tape  data  are  as  follows: 

a.  NGAGE  -  22 

b.  MAXD2  -  5957 
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Typical  compiler  execution  times  (using  optimization  for  fast  execution 
of  object  code,  0PT=2)  and  the  SCM  and  LCM  required  to  load  the  codes  are  given 
in  Table  7. 


TABLE  7.  MEMORY  REQUIREMENTS,  LOAD  OPTION  4 


CODE 

aC0MPILE  TIME 

SCM 

bLCM 

(S) 

(octal  words) 

(octal  words) 

CDEPR0P 

8.383 

132,167 

1,236,310 

DEPROB 

8.335 

111,762 

1,016,020 

a0PT  =  2,  optimization  for  fast  execution  of  object  code. 

^Contains  CBLANK  and  C0M1  and  COM2  common  blocks  with  P  and  Q  set  at  current 
1 imits  (22  x  5957). 

CFor  elastic  analyses  only  (SIGMA  and  SIGMAB  routines  omitted),  SCM  reduces 
to  130,513  words  and  LCM  to  1,132,  360  words. 

Once  the  code  has  been  executed,  an  absolute  object  code  file  is  created 
by  the  segmentation  loader.  Saving  this  file  will  eliminate  the  need  to 
recreate  an  absolute  file  each  time  the  code  is  used,  unless  changes  are  made 
to  either  the  segmentation  directives  or  the  code.  It  is  not  necessary  to 
preset  the  core  to  zero  prior  to  execution. 

Table  8  contains  the  Job  Control  Language  (JCL)  required  to  run  N0VA-2LT 
(N0VA-2LTS)  on  the  CYBER  176  computer,  using  Load  Option  4. 

5.  INPUT  DATA  INSTRUCTIONS 

Most  of  the  required  input  data  pertain  to  DEPROB  or  DEPROP,  and  are 
fully  described  in  References  1  and  6.  The  general  NOVA  input  data  have  been 
altered  and  are  listed  in  Table  9. 

All  the  data  are  organized  in  groups,  with  each  group  beginning  on  a 
separate  data  card.  Additional  cards  may  be  required  for  a  particular  group. 

The  format  corresponding  to  each  group  is  given  in  parentheses  and  is  always 
in  fields  of  12.  Each  parameter  in  the  data  set  is  described  and  its  units 
specified;  the  associated  variable  name  is  also  given.  Columns  73  and  80  can 
be  used  to  label  the  cards  in  the  data  sets  to  facilitate  assembly  of  the 
card  deck  and  recognition  of  the  variables  in  the  data  sets. 

24 


AFWL-TR-79-160 


To  accommodate  the  changes  in  DEPROP  for  the  N0VA2S  version  (Ref.  6), 
the  following  modifications  were  made  in  the  DEPROP  input  data. 

In  Group  5,  page  114  of  Reference  6,  the  variable  NBND  (boundary  condition 
code)  is  defined  as 
NBND  =  PQ 

where  P  and  Q  indicate  the  numerical  codes  for  the  boundary  condition  in  the 
Y  and  8  directions  as  follows: 

Code.  Edge  Condition 

1  clamped-clamped 

2  simple- simple 

3  free- free 

4  clamped- simple 

5  clamped- free 

Values  of  PQ  may  range  from  11  to  55. 

In  Group  17,  page  18  of  Reference  6,  new  groups  A  and  B  are  added  as 
fol lows: 

Group  17A:  (112)  NSPR 

Number  of  linear  elastic  springs.  (NSPR) 

If  NSPR  =  0,  skip  Group  17B. 

Group  17B:  (3112,  F12.1)  IDIR(I),  NSPG(I),  NSPB(I),  BIGK(I) 

Code  designating  direction  of  spring  (IDIR): 

1  u-di recti  on 

2  v-di recti  on 

3  w-direction 

Number  indicating  the  gamma  (Yj)  spring  position.  (NSPG) 

Number  indicating  the  beta  ((3^)  spring  position.  (NSPB) 

Linear  elastic  spring  stiffness  (K,  lb/in).  (BIGK) 

Repeat  Group  17B  for  I  -  1,  NSPR.  The  cards  in  Group  17B  may  be  arranged 
in  any  order. 

Group  23  establishes  a  correspondence  between  the  location  of  the  pressure 
data  channels  on  the  data  tape  and  the  mass  points  of  the  analytical  model. 

The  numbers  which  specify  the  location  of  the  gages  relative  to  the  m^ss  points 
are  ordered  from  the  first  mass  point  for  beam  models  and  from  the  coordinate 

origin  for  panel  models.  As  an  example,  Figure  5  shows  a  simple  beam  model 
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TA8LE  8.  JOB  CONTROL  LANGUAGE 


Job  Control  Sequence 

1  _ 


Comments 


Job  Card 

Account  Card 

ATTACH,  OLDPL, •••  • 

UPDATE  old  program  library  file 

UPDATE,  W,  Q. 

Sequential  library,  quick  compile 

RETURN,  OLDPL. 

FTN ,  A,  I,  LCM=I ,  0PT=2 ,  PL=XXX,  .  .  . 

If  LCM  exceeds  131,071  words, 
LCM=I  is  required.  If  line  print 
limit  exceeds  default,  LP=XXX  is 
requi red. 

RETURN,  COMPILE. 

REQUEST,  A8S,*PF. 

Ensures  proper  file  residence 
for  permanent  file  ABS 

SEGLOAD,  B=ABS. 

LDSET ,  PRESET=ZERO ,  MAP=8SEX,..  . 

LOAD,  LGO. 

NGO. 

REWIND,  ABS. 

REQUEST,  TAPEIO ,  .... 

Pressure  Data  Tape 

ABS. 

Execute  Code 

CATALOG,  ABS,  ...  , 

Absolute  object  code  permanent 
file  for  future  jobs 

EXIT,  S. 

EXIT. 

7/8/9 

UPDATE  Directives 

7/8/9 

Segmentation  Directives 

7/8/9 

Run  Data 

6/7/8/9 


TABLE  9.  N0VA-2LT  (N0VA-2LTS)  INPUT  DATA 
Group  1:  (112)  NCASES 

Number  of  cases  to  be  run. 

Group  2:  (20A4)  TITLE 

Identifying  title 

Group  3:  (5112)  ICOMP,  KTVPE,  KDAM,  KDS ,  NDBUG 

Structural  Element  Position  Code  (ICOMP) 

1  Nonfuselage  element,  or  at  least  no  additonal  skin  effect. 

2  Fuselage  element.  A  stringer,  longeron  or  frame  will 

derive  additional  skin  support  against  internal  pressuriza¬ 
tion.  (See  Vol.  I  of  Ref.  1,  section  4.1.8,  for 
additional  discussion). 

Structural  Element  Code  (KTYPE) 

1  Single  layer  metal  panel 

2  Single  Layer  plastic  panel 

3  Honeycomb  metal  panel 

4  Honeycomb  plastic  panel 

5  Multilayer  plastic  panel 

6  Stringer  or  longeron 

7  Frame 

8  Metal  ring 

9  Plastic  ring  (radome) 

10  Metal  rib 

Damage  Criteria  Code  (KDAM) 

0  No  permanent  damage  (threshold) 

1  Catastrophic  damage 

2  No  criteria,  response  run  only 

100  Noniteration,  but  program  accumulates  maximum  values 

for  threshold  damage  criteria. 

101  Noniteration,  but  program  accumulates  maximum  values 

for  catastrophic  damage  criteria. 

(Note--KDAM  must  be  0,  100,  or  101  if  NL0AD  =  4  in  Group  6). 

Response  option  code  (KDS) 

1  Static  only 

2  Dynamic  only 

3  Static  and  dynamic 

(Note--0nly  the  dynamic  option  is  available  for  KTYPE=10). 
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TABLE  9.  CONTINUED. 


Debug  option  (NDBUG) 

0  No  debug  output  (default  option) 

1  Most  debug  output 

2  All  debug  output 

If  KDAM=2  or  KDS=1,  skip  Group  4 
Group  4:  (F12.1)  ROAM 

Probability  of  exceeding  specified  damage  level,  expressed  as  a 
fraction.  (0  <P0AM  <  1.0) 

For  KTYPE  <  6,  use  DEPROP  input  here.  (Groups  1-14,  Ref.  1) 

For  KTYPE  >  5,  use  DEPROB  input  here.  (Groups  1-22,  Ref.  1) 

If  KPS  =  2,  skip  Group  5. 

Group  5:  (F12.1)  PS 


Uniform  static  preload  pressure,  psi. 

If  KPS  =  1,  skip  Groups  6~20 
Group  6:  (112)  NLOAO 

Dynamic  load  option  (NIOAD) 

1  Uniform  spatially,  linear-exponential  temporal  function. 

2  Uniform  spatially,  point-by  point  temporal  description. 

3  General,  discrete  load  description  (Not  appropriate  for 

KTYPE  =  10). 

4  Data  tape.  Evenly  spaced  times,  one-dimensional  variation 

spatial ly. 


If  NLOAO  =  2,  skip  to  Group  8. 

If  NLOAO  =  3,  skip  to  Group  10  or  15,  depending  on  KTYPE 

If  NLOAO  =  4,  skip  to  Group  21 

Group  7:  (6F12.1)  PP1,  PPO,  TTO,  TPRIME,  AA,  ANN 

Pressure,  Pp  psi  (PP1) 

Pressure,  P  ,  psi  (PPO) 


Time,  t  ,  sec 
Time,  t1 ,  sec 


(TTO) 

(TPRIME) 
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TABLE  9.  CONTINUED. 

Constant,  a,  dimensionless  (AA) 

Constant,  n,  dimensionless  (ANN) 

(Note--See  Figure  1  for  definition  of  parameters). 

Skip  Groups  8-23 
Group  8:  (112)  NTIME 

Number  of  points  to  be  specified  in  point  by  point  load  description. 

Be  sure  to  include  time  =  0.  If  time  exceeds  last  value  in  table,  the 
the  last  value  of  pressure  is  used.  Dimensioned  for  20  times  and 
pressures. 

Group  9:  (2F12.1)  TT(I),  PT(I) 

Time  table,  sec  (TT ( I ) ) 

Pressure  table,  psi  (PT(I)) 

Repeat  Group  9  for  1=1,  NTIME 

Skip  Groups  10-23 

If  KTYPE  <  6,  skip  to  Group  15 

Group  10:  (112)  NPS 

Number  of  pressure  stations  for  which  there  are  valid  data.  Must  have 
at  least  2  stations,  but  not  more  than  41. 

Group  11:  (6F12.1)  SP(I),  I  =  1,  NPS 

Position  of  pressure  station  relative  to  mass  point  locations  in  DEPROB. 
For  example,  a  station  2/3  of  the  distance  between  mass  numbers  3  and  4 
in  the  OEPROB  model  would  mean 

SP  =  3.67.  (0  1  SP  IN  +  1). 

(Note--Do  not  assign  an  SP(I)  to  bad  channels). 

Group  12:  (6112)  LTIME(I) ,  I  =  1,  NPS 

Number  of  entries  in  pressure  versus  time  table  for  each  spatial 
location.  LTIME  must  be  at  least  2,  but  not  more  than  14. 

Group  13:  (6F12.1)  TTB(K,I),  K  =  1,  LTIME(L) 

Break  point  times  for  each  spatial  station,  sec. 
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TABLE  9.  CONTINUED. 

Group  14:  (6F12.1)  PRTB  (K,  I,  LTIME(I) 

Break  point  overpressures  corresponding  to  times  in  Group  13  for  each 
spatial  station,  psi. 

Repeat  Groups  13-14  for  each  station,  1=1,  NPS 
Skip  Groups  15-23 
G  nup  15:  (2112)  NPX,  NPY 

Number  of  pressure  stations  in  gamma  direction.  Must  have  at  least  2 
stations,  but  not  more  than  10-  (NPX) 

Number  of  pressure  stations  in  beta  direction.  Must  have  at  least  2 
stations,  but  not  more  than  10.  (NPY) 

(Note--Since  the  pressure  mesh  is  rectangular,  NPX-NPY  pressure  records 
must  be  supplied). 

Group  16:  (6F12.1)  XP(I),  I  =  1,  NPX 

x-positions  at  which  time  histories  are  specified,  in.  (XP) 

Group  17:  (6F12.1)  YP(I),  I  =  1,  NPY 

y-positions  at  which  time  histories  are  specified,  in  or  deg.  (YP) 

Group  18:  (6F12.1)  KTIME  (J,  I),  J»l,  NPY 

Number  of  entries  in  pressure  versus  time  table  for  each  spatial  point. 
REPEAT  Group  18  for  I  =  1,  NPX 

Group  19:  (6F12.1)  TTP  (K,  J,  I),  K=l,  KTIME  (J,I) 

Break  point  times  for  each  spatial  station,  sec. 

Group  20:  (6F12.1)  PRT(K,J,I),  K=l,  KTIME(J,I) 

Break  point  overpressure  corresponding  to  times  in  Group  19  for  each 
spatial  station,  psi. 

REPEAT  Groups  19  and  20  for  each  station,  with  the  y-index  varyinq  more  rapidl 
J=l,  NPY;  then  1=1,  NPX. 

Skip  Groups  21-23 

Group  21:  (2F12.1)  TIM1,  SKIP 

Start  time,  relative  to  the  data  tape  time  scale,  sec.  Time  should  be  as 
near  as  possible  to  the  first  shock  engulfment. 
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TABLE  9.  CONTINUED. 

Skip  frequency,  dimensionless.  The  time  history  will  be  sampled  at  this 
frequency.  A  1.0  will  mean  every  time  is  used;  3.0  means  every  3rd  time, 
etc.  Must  be  a  whole  number. 

Group  22:  (112)  NGAGE 

Number  of  channels  of  data  on  the  tape.  Each  channel  represents  a 
pressure  gage  and  a  tape  file. 

Group  23:  (6112)  NORDER(I) ,  1=1,  NGAGE 

Gage  location  code.  The  numbers  which  specify  the  location  of  the  gages 
relative  to  the  mass  points  in  the  model  are  ordered  from  beginning  of 
DEPROB  model  (mass  point  1)  or  coordinate  origin  in  DEPROP.  A  zero 
indicates  data  for  a  specific  channel  is  not  to  be  used.  The  example 
in  Figure  5  defines  the  code  for  a  beam  where  the  data  tape  was  generated 
backward  with  four  gages.  Therefore,  the  pressure  data  from  gage  number 
4  is  the  first  data  channel  written  on  the  tape,  and  is  designated  NORDER(l). 
Note  that  gage  number  2  wais  omitted  since  it  was  found  to  be  faulty. 

Table  10  is  constructed  in  the  same  order  as  the  tape  was  written;  the  Ith 
index  of  NORDER  corresponds  to  the  Jth  index  of  SP. 

If  KTYPE  <  5,  the  next  input  is  identical  to  that  described  in  Groups  15,  16, 

and  17.  The  only  limitation  is  that  either  NPX  or  NPY  (or  both)  must  be  equal 

to  1.  If  they  are  both  1,  for  a  uniform  distribution,  then  Groups  16  and  17 
can  be  two  blank  cards. 

If  KTYPE  >  5,  the  next  input  is  identical  to  that  described  in  Groups  10  and  11. 

The  only  difference  is  that  NPS  can  be  1  for  a  uniform  distribution.  In  that 

case  Group  11  can  be  a  blank  card. 

REPEAT  Groups  2-23  for  each  additional  case,  as  specified  in  Group  1.  If  a 
data  tape  was  used  (NL0AD=4),  the  tape  will  be  rewound  and  used  with  the  next 
case  for  which  NL0AD=4  is  specified. 
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with  pressure  gage  locations  and  Table  10  lists  the  pressure  gage  order  and 
locations  relative  to  the  mass  points.  In  this  case  the  data  type  was 
generated  backwards  using  test  data  from  four  gages. 

6.  EXPERIMENTAL  PRESSURE  DATA  (LOGICAL  FILE  TAPE10). 

The  following  procedures  are  suggested  for  generating  an  experimental 
pressure  data  tape  (Logical  File  TAPE10)  for  use  with  N0VA-2LT  (NOVA-2LTS ) : 

a.  Determine,  as  accurately  as  possible,  the  shock  arrival  time  for 
the  first  gage  to  be  engulfed. 

b.  Digitize  data  for  all  gages  approximately  10  us  prior  to  and 

140  ps  following  the  shock  arrival  time  determined  above 

c.  Digitize  all  data  at  10  ps  intervals. 

d.  If  possible,  order  data  on  tape  for  all  gages  used,  assuming 

this  is  consistent  with  the  DEPROB  model.  Not  imperative, 
however. 

e.  Output  data  tape  should  have  time  in  microseconds  ar d  over¬ 

pressure  in  psi.  Tape  should  be  unlabelled,  unformatted, 
stranger(S),  at  800  bpi  using  Block  Type  K  and  Record  Type  U. 

f.  Check  pressure  plots  for  any  suspect  type  gages.  If  data  is  bad 

the  gage  can  be  eliminated  during  execution  of  the  N0VA-2LT 
(N0VA-2LTS)  code  by  means  of  card  input  data. 

g.  Determine  At  for  N0VA-2LT.  Then  select  a  sampling  frequency 

(SKIP)  for  the  data  tape.  For  example,  if  A  =  30  ps  and  the 
tape  is  digitized  every  10  ps,  a  SKIP  parameter  of  at  least 
3.0  should  be  used.  Larger  values  should  be  tried  for  one  case 
to  determine  the  largest  appropriate  SKIP. 

The  format  for  writing  the  digitized  data  on  tape  is  given  in  Table  11. 
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NOTE:  MASS  POINTS  (•)  ARE  ORDERED  COUNTER¬ 
CLOCKWISE  AND  PRESSURE  GAGES  (■) 
CLOCKWISE 


Figure  5.  DEPROB  analytical  model  with  pressure  gage  locations. 


TABLE  10.  PRESSURE  GAGE  ORDER  AND  LOCATION 


Gage  No. 

I=1,NGAGE 

Gage  Order 

NORDER  (I) 

Gage  Location 
aJ  bSP(J) 

1 

3 

1 

C.O 

2 

0 

3 

2 

2 

1.5 

4 

1 

3 

6.0 

aThe  J  indices  (with  the  corresponding  SP‘s)  must  be  ordered  in  increasing 
value. 


^Number  of  gage  locations  (SP's)  must  equal  the  number  of  nonzero  values  in 
the  NOROER(I)  column. 


33 


AFWL-TR-79-160 


TABLE  11.  DATA  TAPE  FORMAT 


Record 

No.  of  Words 

Description  of  Parameter 

EOF 

1 

Identification 

100 

Mixture  of  both  fixed  and  real 
numbers.  Need  to  equivalence 
both  types  of  variables,  e.  g.  , 
ID(100) ,  DD(100) . 

DATA  for 

Channel  1 

ID(17)*ID(18) 

(usual ly 

2*252) 

252  data  pairs  consisting  of  time 
(ps)  and  pressure  (psi),  one  after 
another.  Real-time  numbers. 

Last  Data 

Record  - 
Channel  1 

Same  as  other 
data  records 

Same.  Last  record  may  be  padded 
with  OCTAL  constant  3  followed 
by  nineteen  4's. 

EOF 

1 

Data  for 

Channel  2 

Same  as  other 
data  records 

Same 

Last  Data 
for  Channel 

NGAGE 

Same  as  other 
data  records 

Same 

EOF 

1 

TRAILER 

Record 

10 

Each  word  is  made  up  of  10  F's 

EOF 

1 

EOF 
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IV.  GEOMETRIC  APPLICATIONS  OF  RADIAL  IMPERFECTIONS 

The  DEPROP  modal  solution  in  the  NOVA-2S  computer  code  assumes  that  the 
analyzed  panel  has  a  partial  right-cylindrical  geometry.  For  many  aircraft 
structural  applications,  the  panels  of  interest  are  smooth  but  arbitrarily 
curved.  The  analysis  of  a  smoothly  curved,  noncyl indrical  panel  can  be  attempted 
with  a  circular  approximation  for  which  the  radial  distance  from  the  circle  to 
the  actual  curve  is  minimized.  The  use  of  this  circular  approximation  may  be 
adequate  for  analyses  in  which  large  deformations  are  encountered  (e.  g. , 
catastrophic  failure  analysis),  but  may  prove  to  be  misleading  in  elastic  or 
threshold-of-permanent  damage  analyses  in  which  the  initial  geometry  plays  a 
more  significant  role  in  the  ensuing  deformation  profile  and  time  history. 

For  these  situations  a  more  accurate  model  of  the  panel  curvature  can  be  derived 
by  placing  initial  radial  imperfections  on  the  approximated  circle  to  conform 
with  the  actual  geometry  of  the  panel. 

These  initial  radial  imperfections  are  applied  by  specifying  initial 
preload  amplitudes  for  the  DEPROP  modes  used  to  approximate  the  radial  deflec¬ 
tion  (w-direction)  of  the  analyzed  panel.  The  proper  summation  of  these  initial 
model  amplitudes  will  allow  the  original  circular  model  to  represent  the  actual 
panel  curvature  more  accurately. 

As  a  first  approximation  the  user  may  be  able  to  determine  the  magnitude 
of  the  modal  coefficients  needed  by  superimposing  the  circular  approximation 
on  the  actual  curve  and  measuring  the  amplitude  of  the  major  modal  imperfec¬ 
tions.  For  example,  in  Figure  6,  the  value  6  is  the  amplitude  of  the  second 
symmetric  mode  which  is  the  apparent  dominant  mode  of  the  actual  shape. 

Because  of  the  sign  convention  in  DEPROP,  the  amplitude  6  would  be  assigned  a 
negative  value  (for  the  second  symmetric  mode)  in  the  input  data;  i.e., 

FG(N,M),  Group  13  in  N0VA-2LT  and  Group  30  in  N0VA-2LTS. 
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Figure  6.  Panel  curvature,  idealized  vs  actual 
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The  individual  modal  amplitudes  may  be  calculated  more  accurately  by  the 
utilization  of  a  Fourier  analysis  as  presented  in  Reference  8.  The  Fourier 
method  applies  a  least-squares  fit  to  both  the  actual  and  the  approximating 
periodic  functions.  For  the  DEPROP  analysis,  the  approximating  periodic 
functions  are  the  assumed  modes  defined  over  the  interval  (0  to  n)  as  given 
in  Reference  1. 

First,  the  approximating  circle  must  be  fitted  to  the  actual  panel 
curvature.  This  may  be  done  by  solving  the  following  simultaneous  equation 

for  V  ,  Z  ,  and  R  in  the  Y-Z  rectilinear  coordinate  system. 

o  o 


(V0  -  Vx)2  -  (ZQ 


ZJ2  =  R2 


(y0  *  y2)2  *  (Z0  *  Z2)2  =  R2 
(Y_  '  Y,)2  +  (Z„  -  Z,)2  =  R2 


where 


Y2 ,  Z2 
Y3 ,  Z3 


Coordinates  of  the  circle's  origin. 

Coordinates  of  the  two  end  points  of  the  actual  curve. 


User  determined  coordinates  of  a  third  point  to  be  located 
on  the  approximating  circle. 

Radius  of  approximating  circle. 


The  approximating  circle  is  now  used  to  determine  the  amplitudes  of  the 
initial  modal  imperfections  through  the  use  of  a  Fourier  analysis.  A  function 
f(p)  is  defined  as  the  radial  distance  from  the  circle  to  the  actual  curve 
(radially  inward  distances  are  positive).  This  function  is  then  approximated 
by  a  summation  of  periodic  functions,  G^(p),  and  their  individual  amplitudes, 


f(P)  =  I  A.  G.  (P)  (6) 
k=l  K  K 

The  values  of  A^  are  determined  by  applying  a  least-squares  fit  over  the 
interval  0  to  tt  (the  interval  over  which  the  DEPROP  modes  are  defined)  in  the 
fol lowing  manner: 
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'V  T  N  l2 

f  f(p)  -  I  A.  G,  (3)  dp  = 

Jo  L  k=l  k  <  j 


irn  mmum 


from  which  the  following  condition  must  be  satisfied: 

N  \  /  N 


/oV(p) '  UAkVP)A*=iGj(p) ) dp  =  0 


When  the  orthogonality  constraint 
n 


/oGJ 


(3)  G  (3)  dp  =  0  for  j  f  k 


(7) 


(8) 


is  introduced  into  Equation  8,  the  following  relation  is  obtained 

7T„  71 . 


I  Gk‘«  f<f> \  J0hw)‘ 


dp 


(9) 


However,  the  function  f (3 )  is  not  easily  defined  over  the  continuous 
domain,  requiring  that  a  discrete  domain  approximation  be  applied  by  defining 
f(8r).  The  function  f(3r)  is  a  set  of  values  of  f(p)  tabulated  at  equal 
intervals  3r-  The  range  of  r  is  from  1  to  P  for  which  P-1  independent  values 
of  f(3)  are  determined.  For  the  end  points  it  is  assumed  that  f(l)  =  f(P)  =  0. 


Equation  10  then  takes  the  form 

P-1  P-1 

1  ,  Gk  (  O^r)  '  \  Z, 

r=l  r=l 


and  each  value  of  A^  is  given  by 

P-1 

G.  (r  \f It 

r 


\  = 


I 

r=l 


G.  (3  if(P 
k\'  r)  V  i 


P-1 
I 
r 


:=1  (Gk{^r))2 


(For  k  =  1  to  N) 


(10) 


where  N  is  the  number  of  modes  used  in  the  approximation. 

The  orthogonal  modes  are  defined  below  for  each  set  of  boundary  conditions 
in  0EPR0P . 

a.  Simply  supported  at  both  ends: 

W  =  sin(kPr) 
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b.  Clamped-clamped 


ckor)  =  cosh(-ri)  - 

where 


a.  sinh 
k 


A.p 

kKr 


:  roots  of  cos  Ak  cosh  Ak  =  1 

cosh  A-k  -  cos  \k 
ak  ~  sinh  Ak  -  sin  Ak 

Using  the  following  trigonometric  identities, 


sinh  z  = 


ez  -  e-z 


cosh  z  =■ 


ez  +  e-z 


-  sin 


kHr 


(ID 


and  defining  _  A. 

~  JT 


W  - 


l^'a k\  ^k^r  /  ^  +  ak\  -Ak^r  —  — 

( - j  e  ( - )  e  -  cos  Aup„  +  au  sin  A^p 


kMr  “k 


kMr 


(12) 


2  '  2 
c.  Clamped- simply  supported 

The  same  as  cl amped- clamped  except  that  the  definition  of  Ak 


is  as  follows: 


Ak  :  roots  of  tan  Ak  =  tanh  Ak 

The  modes  are  all  independent  of  each  other  so  that  the  calculation  of  Ak 
for  one  value  of  k  is  independent  of  any  other  value  of  k.  Because  of  the 
characteristics  of  the  least-squares  fit  and  to  the  numerical  accuracies  of 
the  iterated  values  for  Ak  and  ak,  a  fine  mesh  of  points  must  be  chosen  for  Pf. 
An  amplitude  should  be  calculated  for  each  mode  used  in  the  analysis,  and  the 
values  of  Ak  back  substituted  into  the  initial  equation. 


N 

fOr)  =  2  Ak  Gk(Pr)  (13) 

to  determine  the  loss  of  accuracy  due  to  numerical  round  offs. 

Note  that  the  calculated  amplitudes  of  the  initial  imperfections,  Ak> 
referred  only  to  the  modes  in  the  £  direction.  The  calculated  values  of  Ak 
are  therefore  associated  with  the  coefficients  of  the  first  A  mode  and  the 

38 


AFWl-TR-79-160 


kth  p  mode.  This  will  give  an  accurate  representation  of  the  actual  panel 
curvature  only  at  the  center  of  the  assumed  model,  and  the  accuracy  will 
decrease  as  the  A  boundaries  are  approached  However,  the  analysis  can  be 
applied  to  the  entire  surface  of  the  panel  and  include  both  the  A  and  p 
imperfection  modes.  The  A  modes  would  be  used  to  preserve  an  approximate 
straight  shape  in  the  A  direction  over  the  noncircular  cross-section,  except 
at  the  A  boundaries  where  by  definition  the  imperfections  are  zero.  The  initial 
approximation,  presented  in  Equation  6,  is  more  generally  given  by  the  following 
equation: 

'Mr)1  "k=1  J.yiNWS)  (14> 

By  applying  the  analogous  solution  procedure  given  in  Equations  7-10,  the  value 
of  the  initial  imperfection  modal  coefficients,  A£k,  are  given  in  Equation  15. 


P-1  Q-l  P-1 

*=1  GkWH»(S)  '(VMNk  ^ 


Q-l 

i 

s=l 


(Mr'"*  (0/ 


(15) 


for  which  H£(A$ ^  is  the  2th  mode  in  the  A  direction,  and  it  is  analogous  to  the 
p  modes  described  above  for  the  applicable  boundary  conditions. 

The  A  function  is  divided  into  equal  intervals  such  that  s  =  1  to  Q  and  A^ 
is  the  coefficient  of  the  2th  A  mode  and  the  kth  p  mode. 
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APPENDIX  A 

DYNAMIC  LOAD  OPTION  CHANGES 

This  appendix  contains  the  changes  necessary  to  convert  the  NOVA-2  code 
to  the  N0VA-2LT  version  and  replaces  the  nuclear  blast  loading  model  with 
four  transient  pressure  models  and  one  static  pressure  model. 

Since  the  NOVA-2  code  library  is  maintained  at  the  AFWL  using  the 
Control  Data  Corporation  6000/7000  UPDATE  program,  the  UPDATE  directives 
needed  to  accomplish  these  changes  are  also  included. 
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CHANGES  tc  CONVERT 
NOVA-2  TO  NCVA-2LT  TAPE  VERSION 


* I DENT  JAN0179 
•DELETE  LIST1.158 

1  11X»5HV(IN)  *11X.5H»  UN)  .7X.14HPRESSURE  (PSD) 

•OELETE  LIST2.165 

1  11X.5HV  (IN)  *  1 1 X  ,  51-  V*  <  IN  >  .7X,14HPRESSURE  (PSD) 

•INSERT  SIGMA. 17 

LEVEL  2,  ALTT 
•PUROECK  IODUM 
•PURGE  IOOUM 
•PUROECK  CLOAO 
•PURGE  CLOAO 
•AODFILE  INPUT. CNOVA 
•COMDECK  CLOAO 

COMMON  /CLOAO/  PP 1 . PPC , TTC , TPP I ME . A A . ANN ♦ OTT 1 . OT TC . A Z . 

1  JL.NTIME.NLOAC.PT (20) ,TT f20) .  KTIME(lO.lO)  LTIME (AD . ISP (40) . 

2  JL8 ( 4 D »NPS . TTP ( 6 . 1 0  » 1 0 ) .oR T ( 6 . 1 0 ♦ 1 0 ) .NPX .NPY , XP ( 22 ) ,YP(22) . 

3  IX I  (23)  . JYJ(23) , JLT (10.10) .PRTT(IO.IO) .DXl (23) .0Y1  (23) , 

4  NGSUM.CEL.MAXC2.MAXC.PS 
•COMCECK  COM1 

COMMON  /CCmi/  P  ( 22 . 5*55  7 ) 

LEVEL  2.  P 
•COMCECK  COM2 

COMMON  /CCm2/  0(22.5957) 

LEVEL  2.  G 
•PURCECK  NOVA 
•PURGE  NOVA 
•AOOEILF  INPUT.CRLANK 
•OECK  NOVA 

PCOGPAM  N C v A  < INFUT»CuTPLT.TAPE5=INPUT,TAPE6=OLTPLT,TAPE1=513. 

1  T APE  1 0=  D 

THIS  IS  THE  N0VA-2LT  VERSION  OF  NOVA.  THE  AERODYNAMIC  AND  BLAST 
ROUTINES  are  PEPL«CEC  eY  USER  DESIGNATED  PRESSURE  FUNCTIONS. 
PROVISION  HAS  ALSO  BEEN  ma-'E  FOR  READING  PRESSURE  CATA  FROM  TAPE. 
NOVEMBER. 1977. 

CALL  CLOAO 
•CALL  CNOVA 
C 

1  FORMAT (6112) 

2  FORMAT (6F 12. D 

3  FORMAT  ( 2 0 A4 ) 

NCASE  =  0 
INOUT  =  1 
REA0(5.1)  NCASES 

100  REA0(5,3)  (TITLED)  .1*1.20) 

NCASE  *  NCASE  ♦  1 
KERR  =  0 
NTRIAL  =  0 

RE  AC  (5.D  ICOMP.KTYFE .KCAm.KDS.NDBUG 
IF  (KDS.EC . 1 )  KC AM  s  2 


42 


NCHPT  S  o 

IF  (KDAM.LT ,?)  NCHPT  =  1 
IF  (KDAM.GT.2)  KC AM  =  KOAM  -  100 
IF (KDAM.LT ,2)  RE AC  (5*2)  PrAM 
IF (INOUT.EO.O)  GC  TC  1400 
WRITE (6.3000)  (TITLE <I) , Is?, 20) 

IF  (KTYPE.LT.6.CR.KTTPE.GT.7)  ICOMP  *  5 
IF  ( ICOMP. EO. 2)  WRITE  (6,5o0O) 

IF  (NCHPT. EG. 1)  WRITE  (6*3100) 

IF  (NiCHPT.EO.fl)  WRITE  <6,3?00) 

IF (KOAM.EQ.O)  WRITE(6.3300)PDAM 
IF (KOAM. EG. 1 )  WRITE (6*34 00) EDAM 

GO  TO  (300, 400. 500. 600,700,800. 900. 1010, lOOn, 1020) *KTYPE 
300  WRITE (6.3500) 

GC  TO  1050 
<*00  WRITE  (6,3600) 

GC  TO  1050 
500  WRITE (6,3700) 

GC  TO  1050 
600  WRITE (6.3800) 

GO  TO  1050 
700  WRITE (6,3900) 

GO  TO  1050 
800  WRITE (6.4000) 

GO  TO  1050 
900  WRITE (6.4100) 

GO  TO  1050 
1000  WRITE (6.4200) 

GO  TO  1050 
1010  WRITE (6,4700) 

GO  TO  1050 
1020  WRITE (6,4800) 

1050  GO  TO  (1100,1200,1300),  KC? 

1100  WRITE (6.4300) 

GO  TO  1400 
1200  WRITE (6,4400) 

GC  TO  1400 
1300  WRITE (6,4500) 

1400  NCAll  =  2 

IF (KTYPE.GT .5)  CALL  CEpRCB 
IF  (KTYPE.LT. 6)  CALL  CEPRCc 
IF (KERR.GT.O)  GC  TO  16qo 
NCALL  =  1 
CALL  PINIT(O) 

IF (KTYPE.GT .5)  CALL  CEPRC? 

IF  (KTYPE.LT. 6)  CALL  CEPRCP 
TF(KDS.EC.l)  GO  TC  1600 
IF (KERR.GT.O)  GC  TO  1600 
NCALL  *  0 
KCK  =  0 
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CALL  P I N I T  <  1 ) 

IF  (KEPR.GT.O)  GO  TC  1700 
RTRIAL ( 1 ) =1 .0 
1500  NTRIAL  =  MPIAL  ♦  1 

IF (KDAM.LT.2)  WRITE (6*4600)  NCASE ,NTRI AL .RTRI AL  ( 1 > 

IF  (KTYPE.GT.5)  Call  CEPPCn 
IF  (KTYPF.LT. 6)  CALL  CEPRCc 
IF (KEPO.NE.0)  GC  TO  1600 
IF  (NCHPT.EO.m  GC  TC  1600 
CALL  PITER  (CRIT ,PTRIAL»NTRIAL*8,K0K) 

IF(KOK.EC.O)  GO  TC  1 5 0 0 
1600  IF  (NCASE. LT.NCASES)  Go  TC  100 
1700  STOP 

c 

3000  FORMAT  ( 1 H 1 ,30X* 15HN  C  V  A  -  2  L  T//1X.20A4) 

3100  FORMAT  ( 1 Ah  ITERATION  PUN) 

3200  FOPMAT  ( 32h  RFSFCNSE  PUN  ONLY*  NO  ITERATION) 

3300  FORMAT  (A2h  NO  DAMAGE  LEVEL*  PROBABILITY  OF  EXCEEDING*  Ffc.3) 

3400  FORMAT  (52h  CATASTROPHIC  DAMAGE  LEVEL*  PROBABILITY  CF  EXCEEDING, 
1F6.3) 

3500  FORMAT  (28h0SINGLE-L^ TER  McT AL  PANEL  ) 

3600  FORMAT  (30H0SINGLE-LAYEP  PLASTIC  PANEL  ) 

3700  FORMAT  (25ho HONEYCOMB  METAi  PANEL  ) 

3800  FORMAT  (27H0H0NEYC0Me  PLASTIC  PANEL  ) 

3900  FORMAT  (29H0MULTI-LAYER  PLASTIC  PANEL  ) 

AOOO  FORMAT  (30HOMETAL  STRINGER  CP  LONGERON  ) 

4100  FORMAT  (15H0MFTAL  FRAME  ) 

4200  FORMAT  (16H0PLASTIC  PING  ) 

4300  FORMAT  (21H0STATIC  SCLUTICn  ONLY) 

4400  FORMAT  (22HOOYNAMTC  RESPONSE  ONLY) 

4500  FORMAT  (37H0STATIC  SCLUTICn  AND  DYNAMIC  RESPONSE) 

4600  FOPMAT  (12H1CASE  NUMBER  12/ 

114H  TRIAL  NUMBER  13,  lOX,  j€H  RANGE  FACTOR  =  E14.6) 

4700  FORMAT  (  UMOMETAL  RING) 

4800  FORMAT ( 13H0RIB  BUCKLING) 

5000  FORMAT  (69H0STRLCTURAL  ELEmENT  DOES  OERIVE  ADDITIONAL  SUPPORT  FROM 
1  FUSELAGE  SKIN) 

END 

•PURDECK  PINIT 
•PURGE  PINIT 
•AOOFILE  INPUT  *  CSETUP 
*OECK  PINIT 

SUBROUTINE  RINIT(M) 

•CALL  CNOVA 
•C*LL  CLOAD 

•CALL— CPlsK  1-  -  ........  —  . .  ......... - - 

•CALL  CCM1 
•CALL  COM2 

DIMENSION  ID  ( 1  00)  *DC  dOO)  *  4  1  ( 1004)  »NORCER  (2?) 

EQUIVALENCE  ( TO  ( l ) »CC  ( 1 ) 


ooo  r>  nr>n  r>or> 


CIMEN5I0N  TTP ( 14,4 1 >  *PRT8< 14,41) »SP (41 > ,SPS  (4 1 ) 
EQUIVALENCE  (PRT<1,1«1),PRtB(1,1)),  (TTP (1,1,1) «  TT8 ( 1 ,  1 ) ) 
DATA  TRO/lOHFFFFFFFFFF/ 

IF(M.EQ.l)  G 0  TC  200 
PS=0.0 

IF (KDS.EC.2)  GO  TC  lSO 
STATIC 

READ (5*2000)  PS 
WRITE (6,2200)  PS 
NU  =  1 
PPP=PS 

IF  (KTYPE.LT. 6)  C*C  TC  ISO 
IF  (KTYPE.LT. 10)  GO  TC  50 
DO  30  I=1*NMASS 
30  PP  (I)  =  0. 

GO  TO  150 

50  DO  100  I=1,NMaSS 
100  PP(I)  s  PS 
150  RETURN 

DYNAMIC 

200  IF(KOS.EC.l)  GO  TC  400 
READ  (5,2050)  NLCAD 
WRITE  (6,2400)  NLCAC 
GO  TO  (250,500,800,6000),  \LOAD 
250  REAO (5,2000)  PP 1 ,PPC , Tto , TpR IME , A A , ANN 
WRITE (6,2300)  PP1 ,PPC , TTC , TPRIME , A A ♦ ANN 
NU=1 

IF (TPRIME. EQ. 0.0)  GC  To  300 
opRIIME=PP0*(1.0  -  TPPlME/TTO)**ANN 
P°R IME  =  PPRIME*EXP (-Aa*TPpIME/TTC> 

TT1=TPRIME*PP1/ (Pc 1-PPRlVE ) 

0TT1=1.0/TT1 
300  OTTOal.O/TTO 

A2=AA*OTTC 
400  RETURN 

500  REAC  (5,2050)  NT  I*E 

REAO (5,2100)  (TT(I) ,PT(I) ,t=1,NTIme) 

WRITE  (6,2500)  NTIME, <TT ( I ) ,PT ( I ) ,1*1, NTIME) 

JL‘  =  Z  .  ' 

NL  =  1 
RETURN 

e00  TF  (KTYPE.GT.5)  GC  TC  1000 
PANELS. 
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RE  AO  (5.2050)  NPX.NPY 

WRITE  (6. 2700)  NPX.NPY 

RE  AO  (5.2000)  < XP (I) ♦ Is  1  ,NP X  ) 

WRITE  (6.3100)  { XF  < I )  .1=1, nFX) 

READ  (5.2000)  <YP ( J) »»*1 ,NRY ) 

WRITE  (6.3200)  (YP ( J > ♦ J*1 .NPY) 

WRITE  (6.3300) 

DC  820  1*1, NPX 

READ  (5,2050)  (KTIME <-. I ) ♦ J=1 .NPY) 

820  WRITE  (6.2800)  ( K  T  I V  E  <J»I)j=l,NPY) 

DC  840  1  =  1  .NPX 
CC  840  J= 1 .MPY 
NTIHE  =  KTIME(J.I) 

RE  AO  (5,2000)  (TTP <K,J.I> ,K  =  1.NTIME) 

WRITE  (6.3600)  I  .^  ♦  ( TTP (K . j , I ) ,K  =  1 .NTIME ) 

WRITE  (6.2900) 

REAO  (5,2000)  (PRT (K »  J,  I ) ,K=1 .NTIME) 

840  WRITE  (6,3000)  (PRT <K « J.I ) ,K=l .NTIME) 

SPATIAL  INTERPOLATION-EXTRAPOLATION.  INDICES  APE  LCWER  8CUN0 
00  900  1=1. NOT 
00  860  III  =  l.NPX 
IE  (XP(III) .GT.XG(I) >  GO  Tr  8«0 
860  CONTINUE 
III  =  NPX 

880  IE  (III.6T.1)  III  *  III  -  1 

0X(I)  *  <  X  G  < I )  -  XP(III) )/(XP(III*l)  -  XP(III)) 

900  I X I (I)  =  III 

CC  960  J  =  l.NRT 
DC  920  JJJ  =  l.NPY 
IF  (YP(JJJ) .GT.XP(J) )  GO  T 0  940 
920  CONTINUE 
JJJ  =  NPY 

940  IE  (JJJ.GT.l)  JJJ  =  JJJ  -1 

0Y1(J)  =  <  XB ( J ) )  -  YP (JJJ) j/(YP(JJJ*l>  -  YP(JJJ)) 

960  JYJ(J)  =  JJJ 
NL  =  0 

DC  980  1=1, NPX 
DC  980  J* l.NPY 
980  JLT(J.I)  =  2 

return 

REAMS. 

1000  READ  (5,2050)  NPS 
WRITE  (6.3400)  NPS 
REAO  (5.2000)  (SPfl) » 1=1 .NPS ) 

WRITE  (6.3500)  ( SF < I ) *  I  =  INnS ) 

WRITE  (6,3300) 

REAO  (5,2050)  (LTIME U) .IslNPS) 

WRITE  (6,2800)  (LT IVE ( I ) » 1= 1 ,NPS) 

CO  1200  1=1, NPS 


NTIME  =  LTIME(I) 


r>  o  o  o 


PE  AO  (5*20  00)  (TTB  (K,  I)  ,K  =  1  ,iKTIME) 

WRITE  (6*3700)  I , ( TTE (K ♦ 1 > = I .NT IME ) 

PESO  (5.2000)  (PRTH(K.I) .Ksl.NTIME) 

WRITE  (6.3800) 

1200  WRITE  (6.3000)  (PRTE (* . I ) .*=1 .NTI^E) 

C  SPATIAL  INTERPOLATION  -  EXTRAPOLATION.  INDICES  ARE  LOWER  80UN0 . 

1 250  DO  1500  1  =  1  . NPS 
ESP  =  SP(I) 

II  *  ESP  ♦  .00001 
Oil  =  ESP  -  FLOAT  (II) 

SPSX  =  0.0 

IE  (II.EC.O)  GO  TC  1^00 
DC  1300  J=1,IT 
1300  SPSX  a  SPSX  ♦  DSOCCJ) 

1400  IE  (II.LT.NMASS*1)  SPS(I)  a  SPSX  ♦  DI  I*OSOC  ( 1 1  ♦  1 ) 

1500  CONTINUE 

00  1600  J=2.NMASS 
1600  CSOO(J)  =  DSOO(J-l)  ♦  DSCC('w) 

DO  1850  1=1 .NMASS 
DO  1700  III  =  1  .NPS 
IF  (SPS(IH)  .GT.DSOCd)  )  C-C  TO  1800 
1700  CONTINUE 

III  *  NPS 

1800  IF  (III.GT.l)  HI  =  III  -  1 

DSOO(I)  =  (OSOO(I)  -  SPS(III) ) ✓ ( SPS  C 1 1 1 ♦ 1 )  -  SPS(III)) 

1850  ISP(I)  =  III 

IF  (MLOAC.EQ.4)  GC  TC  7100 
CD  1  POO  1  =  1. NPS 
1900  JLP ( I )  =2 
RETURN 

LC AO  OPTION  4  -  TAPE  INPUT  FROM  TAPE10. 

PROGRAMMED  FOP  7600  CNLY. 

6000  MflXD  =  11914 
M AX02  =  59S7 
RE AO  (5.2000)  TIM1.SKIP 
T IM2  =  TIM  ♦  TSTLP  ♦  DELTtY 
NSKIP  =  SKIP  «.  .0001 
READ  (5,2050)  NG AGE 
READ  (5,2050)  < NCPDEP < I ) , I = 1 . NGAGE ) 

WRITE  (6,3900)  TIM, NSKIP  .nGAGE  »  (NOROER  ( I )  ♦  1*1  .NGAGE ) 

NGSUM  =  0 
CO  6050  I  =  . NGAGE 

IF  (NOROER(I)  .GT.C)  N G SUM  =  NGSUM  *1 
6050  CONTINUE 

IF  (KTYPE.GT.5)  GC  TC  7000 
C  PANELS. 

READ  (5.2050)  NPX.NPY 
WRITE  (6,2700)  NPX.NPY 


IF  (NPX.GT.1.4ND.KPV.GT.1 )  GO  TO  8800 
BF  AO  (5, 2000)  ( XP ( I ) *  I  =  1 ,NP X ) 

BEAD  <5*2000)  (YP(I) iI=1*NpY) 

V*R I TE  (6*2100)  <XP(J)  *1  =  1, kPX) 

Vk  P I  TE  (6*2200)  (YF(I)  »I  =  1,'FY) 

JL  =2 
MJ  =  1 

IF  (NPX*NPY.EO. 1 )  GC  To  71n0 
Nl'  =  0 

IF  (NPX*NPY .  NE . N6eUM )  GO  Tc  8600 
IF  (NPX.EQ.l)  GC  TO  6600 
CO  6500  1*1, NOT 
CO  6300  1 1 1  =  1  , NPX 
IF  (XP(III) .GT.XG(I) )  GO  To  6400 
6200  CONTINUE 
III  =  NPX 

6400  IF  (III.GT.l)  HI  =  III  -  l 

0X1(1)  =  (XG(I)  -  XP(IlI) ) /(XP(III*1)  -  XP(III)) 
6500  I  X I  ( I )  =  in 
6600  IF  (NPY.EC.l)  GC  TO  7100 
00  6900  J=1,NRT 

CO  6700  JJJ  =  1 , NF Y 
IF  (YP< JJJ) .GT.XB ( J) )  GO  TC  6800 
6700  CONTINUE 
JJJ  =  NPY 

6800  IF  (JJJ.C-T.l)  JJJ  s  JvJ  -  i 

Ol'Yl(J)  =  (  X8  <  J )  -  YP  (JJJ)  )/ (YP(JJJ*1)  -  YP  ( JJJ )  ) 
6900  JYJ(J)  =  JJJ 
GO  TO  7 1 C  0 
C  BEAMS, 

7000  BEAO  (5,2050)  NPS 
ViB  I  TE  (6*2400)  NPS 
BEAO  (5*2000)  (SP(I)  » 1  =  1  *NpS) 

'NRITE  (6*2500)  (SP  ( I  >  *  1  =  1 *nPS) 

JL  =  2 
MJ  =  1 

IF  (NPS. EC. 1)  GC  TO  7 1 00 
IP  (NPS.NE.NGSUM)  GC  To  05nO 
MJ  =  0 
GC  TO  1250 
C 

7100  CEL  *  0 

TCV  =  l.E-6 

- •  -TTMl.  «-TIKi./T«V.  ...  «. . .  .  .  , 

kg  =  0 

KK  =  0 
KKK  s  o 
NTIME  =  0 

BLPFEP  IN  (10*1)  ( IC ( i ) * IC ( 100) 

IF  (UN  I T  (  1  0 ) )  7200*72  00,P3oC 


7200  PUFFER  IN  ( 1 0 « 1 )  ( IC U ) * IC ( 1 00 ) ) 

IF  (UNIT ( 10) )  725C*8lOo»83nQ 
7250  IF  (ID(l).EO.TPC)  GC  To  87r0 
NWORDS  =  10(17) 

NPOINT  s  10(18) 

LR  *  NWOPCS*NPOINT 
NWS  =  NWC«OS*N5KIP 
KK  =  KK  ♦  1 
KG  =  NOPCER(KK) 

IF  (KG.GT.O)  KKK  =  KKK  ♦  l 
IL  a  NWOPOS 

7300  PUFFER  IN  (10.1)  (AI (1) ,AI (LR) ) 

IF  (UNIT  (10))  7400.fi  lOO, 8-5C0 
7400  IF  (DEL.fiC.o.)  CEL  a  < A I ( 1 DWORDS )  -  A I ( 1 ) ) *SP IP*TC V 
C  LOCATE  FIRST  TIME. 

00  7500  IaTL.LR.NWCRCS 
IF  (AI ( I - 1 ) .GE.TIV1)  GO  TO  7600 
7500  CONTINUE 

GO  TO  7300 

7600  IF  (KG.GT.O)  P(KG,1)  =  JIlj) 

T 1  =  AI ( 1-1 ) *TCV 

IF  (NITME.EQ.O)  NTIME  =  ( T  I’M  2-T IM  1  *TCV ) /DEL  ♦  2 
IF  (NTIME. GT.MAXC)  CC  TO  8400 

J  a  1 

IF  (KG.EC.O)  GO  TO  7fi00 
IL  a  NWS  -  LR  *  I 
IF  (IL.GT.O)  GO  TC  7600 
IL  =  I  ♦  NWS 
OC  7700  IslL.LP.NWS 
IX  =  I 
J  =  J  ♦  1 

IF  (J.GT.NITME)  GC  TC  7700 
P(KG,J)  a  AI(I) 

7700  CONTINUE 

7750  IL  =  NWS  -  LR  ♦  IX 

7800  PUFFER  IN  (10,1)  ( A  I ( 1 ) , A  I  f LR ) ) 

IF  (UNIT (10))  7900,8100, 8300 
7900  IF  (KG.EC.O)  GO  TC  7fiC0 

IF  (J.GE.NTIMF)  GC  TC  7800 

DC  8000  I*TL,LR,NWS 
IX  a  I 
J  =  J  ♦  1 

IF  (J.GT.NTIMP*  GC  TC  8000 

IF  (J.GT.MAX02)  GC  TC  7950 

P (KG.J)  =  AI(I) 

GO  TO  8000 

7950  C  (KG.J-MAX02)  =  AKI) 

8000  CONTINUE 

GO  TO  7750 
ENO  OF  GAC-E  DATA. 
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C 


o  o  r> 


fllOO  IF  (NO0UG.GT.O)  WRITS  (6.4?C0)  KK , KG, DEL ♦ T 1 ♦ NTIWE ♦ 
1  NWORDS.NPOINT, (IC (I) »I=1,P> 

IF  (KKK.lT.NGSUN)  GC  Tq  72oC 
C  DATA  REAC. 

8200  REWIND  10 

WRITE  (6.4800)  CEl.NTiME 

IF  (NDRUG.LT. 2)  GC  TC  8250 

DC  8220  I=1.MGSLV 

NT2  =  NTI*E  -  maXC2 

WRITE  (6.4900)  I.  (P (I.J) J*1 .NT1) 

IF  (NT2.GT.O)  WRITE  (6,4800)  I , (Q ( I . J) ♦ J=1 .NT2) 
8220  CONTINUE 
8250  RETURN 

ERRORS. 


8300  WRITE  (6.4300) 
GC  TO  8900 
8400  WRITE  (6.4400) 
GC  TO  8900 
8500  WRITE  (6,4500) 
GC  TO  8900 
8600  WRITE  (6,4500) 
GC  TO  8900 
8700  WRITE  (6.4700) 
GO  TO  8900 
8800  WRITE  (6,4000) 
C 

8900  KERR  a  2 
RFTURN 


ntime.kaxc 

NPS .NGSum 
N P X  , N P T  tNGSl.W 
KK .NGAGe 


2000  FCPMAT(6F12.1) 

2050  FORMAT  (6112) 

2100  FORMAT  2F12.1) 

2200  FORMAT  (24HOSTATIC  PRESSURE,  dsi  =  E15.6) 
2300  FORMAT  (23HODYNAWIC  LCaq  CONSTANTS/ 


11H 
11H 
1  1 M 
11H 
1  1 M 
1  1 M 


PP1 

PPC 

TTO 

TPR I ME 
A  A 
ANN 


El5.fi/ 

El5.fi/ 

El5.fi/ 

El5.fi/ 

El5.fi/ 

El5.fi) 


2400  FORMAT  (21HOOYNAMIC  ICaD  CpTION  14) 
2500  FORMAT  ( 1 8HONUMBER  CE  TIMES  =  I4/28H 
1  (2E15.fi)) 

2700  FORMAT  (24HONUM8ER  CE  LOAD  STATIONS/ 
1  12H  NPX  =  1 3/ 1 2H  Ncy  =  13) 


TIME*  SEC 


2800  FORMAT 
2900  FORMAT 
3000  FORMAT 
3100  FORMAT 


(5X.10I5) 

( 1 8h  PRESSURES  (PSI)  =) 

(5X.6E15.fi) 

( 1 9 HO X- RCS I T I CN  s  (In)  =/ (5X ,5E 15 . 6 ) ) 


PRESSURE,  PSI/ 


3200  FORMAT  (26HOY-PCSITICNS  <IN>  OR  OFG)  =  /SX , 5E 15 . 6  )  ) 

3300  FORMAT  ( 26HONUMBER  CF  TABLr  ENTRIES  =) 

3400  F ORM A T  (24HONUMBER  CF  LOAC  STATIGNS/12H  NPS  =  13) 

3500  FORMAT  ( 25HOMF  ASUREmEN  T  POSITIONS  =  /5X.5E15.6)) 

3600  FORMAT  (14H0TIMES  (SEC)  ,l0X,6HNPX  =  2C.5X6HNPY  a  13/ 

1  <5X,6El5.6)) 

3700  FORMAT  (14H0TIMES  (SEC)  , 1  OX ,2H I  =  1 3/ ( SX ,6E 1 5 . 6 ) ) 

3600  FORMAT  (1<5h  ppeSSLRES  (SEC)  =  ) 

3ROO  FORMAT  (SHOTAPE  USE/  34H  STAPT  TIME.  SEC  (TImD  =  E15.6/ 

1  34H  SKIP  FRECUENCy  (NSKIP)  =  16/ 

2  34M  NC,  OF  GAGES  Cn  TAPE  ( NGAGE )  =  16/ 

3  2SH  LCCaTION  IC  CF  GAGEc  =  /(5X.10I4)) 

4000  FORMAT  (33HOTAPE  INPUT  IS  *NLY  1  DIMENSICNAU/ 

1  28H  EITHEP  NPX  CR  NPy  MUST  BE  1) 

4200  FORMAT  (1SH00ATA  FOR  CAGE  vC ,  14,  15H,  LOCATION  IC  14/ 

1  1 4H  TIME  INTERVAL  ElS.6,  13H.  START  T I ME  E15.6/ 

2  16H  NUMBER  CF  TTMES  16/  ?6h  NUMBER  OF  WOROS  PER  PCINT  14/ 

3  28K  NUMBER  OF  PCINST  PER  RECORD  15/  1X.9A10) 

4300  FORMAT  (26HOPARITY  ERROR  CN  DAT  TAPE) 

4400  FORMAT  (25HOOATA  EXCEEDS  TlBLE  SPACE  215) 

4500  FCRMAT  (32HONUMBER  CF  ACTIyE  GAGES  IS  WRONG  314) 

4700  FORMAT  (12HO  FNC  CF  Tape  214) 

1  1  PH  TIME  INTERVAL  =  E15.6/18H  NO.  CF  TIMES  =  16) 

4P00  FCRMAT  (5MO  I  =  I3.4X.gHP  (PSI)  =  / ( 1 X , 1 OE 1 2 .4 ) ) 

END 

•PUROECK  PRESS 
•PURGE  PRESS 
•AODFTLF  INPUT, INTI 
•CECK  PRESS 

SUBROUTINE  PRESS 
•CALL  CNOVA 
•CALL  CLOAO 
•CALL  CPLK1 
•CALL  COM1 
•CALL  CCM2 

DIMENSION  PPT TP (41) 

DIMENSION  TTP(14.41) »FRTP(14,41) 

EQUIVALENCE  ( pRT (1,1*1)  .PRt6  (1,1)  )  ,  ( TTP (1,1,1), TTE (1,1)) 
EQUIVALENCE  (PRTT ( 1 , 1 ) .PRTtb  ( i ) ) 

C 

IF (NCALL.GT.O)  GC  TC  SOOO 
ZZ  =  1.0/RTRIAL  (1) 

GO  TO  (50.220,800,1000),  NlCaD 
50  IF (TIME.GE.TPRIME)  TC  TO  10* 

PPP=ZZ*PF1*(1.0  -  T IME *OTT 1 ) 

IF (PPP.LT  .0,*)  PPPsO.O 
GC  TO  400 

100  IF (TIME.GE.TTO)  GC  TC  200 

ppo=pPO* (1,0  -  TIME*CTTO) *«ANN 
PPPsZZ*PPR*EXR ( -A Z*T I  ME) 
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GO  TO  400 
200  P PPaO.O 
C-0  TO  400 
C 

220  CC  240  J=JL.NTIMF 

IF  (TIME.LE.TT ( J) )  GC  TO  260 
240  CONTINUE 
JL  =  NTIC'E 
POP  =  7ZapT ( JL ) 

GO  TO  400 
260  JL  =  J 

PPP  =  PT(J-l)  ♦  f  T  I H  E  -  TT  t’v-1 )  )  *  (PT  ( J)  -  PT(J-l))/ 

1  (TT ( J)  -  TT ( J-l ) ) 

PPP  =  22*P°P 
C 

400  IF  (KTYPE.LT. 6)  GC  TC  9000 
PX  a  Poo 

IF  (KTYPE.EQ.10)  FX  =  0. 

CC  500  I=l.NMflSS 
500  PP  (  I)  =  PX 
GC  TO  9000 
C 

POO  IF  (KTYPE.GT.5)  GC  TC  900 
C  PANELS. 

CC  860  1=1, NPX 
CC  860  jsl.NPY 
c  interpolate  ON  TINE. 

PPP  =  0.0 

IF  (TIME.LT.TTP(1,J,I) )  GC  TO  860 
JL  =  JLT(J.I) 

NTIMF  =  KTIME(J.I) 

DC  820  K=JL,NTIME 
KK  =  K 

IF  (TIME.LE.TTP(K,J,I>  >  C-C  TO  840 
820  CONTINUE 

JLT(J.T)  =  NTTve 
PPPsPRT(NTI me , J , I ) 

840  JL  =  KK 

PI  =  PRT(JL-1.J,I) 

T 1  =  TTPUL-ltJ.I) 

PPP  :  PI  ♦  (TIME  -  Tl)*(FRT(JL.J*I)  -  P 1 ) / ( TTP PJL t  w  *  I )  -  T 1 ) 
JLT(J.I)  =  JL 
860  PPTT(J.I)  =  PPP 
C  INTERPOLATE  SPATIALLY. 

K  =  0 

CO  880  1=1, NGT 
II  =  IXI  (I) 

DX  =  DX1  ( I) 

OC  880  J=1»NPT 

IF  <NUSE(J*I) .FC.O)  GC  TC  «80 
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o  o  r> 


K  =  K  ♦  1 
JJ  =  JYJ(J) 

DY  =  0Y1 ( J) 

PI  =  PRTT(jj.II)  ♦  CV*(PRTT(JJ*UII)  -  PRTT(JJ,II)> 

P2  =  PPTT ( J  J » 1 1  ♦  1 )  *  Cy*(PpTT (JJ ♦1.11*1)  -  PRTT ( JJ  *  1 1 ♦ 1 ) ) 
PPP  =  PI  ♦  DX»(P2  -  PI) 

PA (K>  =  PPP»Z7 
880  CONTINUE 
GO  TO  9000 
C  REAMS. 

900  CO  930  1=1. NPS 
PPP  =  0.0 

IP  (TIME.LT.TT9 ( 1 , I ) >  GO  Tr  930 
JL  =  JL9<I) 

NTIME  =  LTIME(I) 

CO  910  K  =  JL. NTIME 
KK  =  K 

IE  (TIME.LE.TTB(K.I) )  GO  TO  920 
910  CONTINUE 

JLP(I)  =  NTI^F 
PPP  s  PRT0 (NTIMF,  I) 

GO  TO  930 
920  JL  =  KK 

PI  =  PRT0UL-1.I) 

T 1  =  TTB(JL-l.I) 

PFP  s  PI  •  (TIME-T1 ) • (PRTe CwL. I)  -  P1/(TT9(JL,I)  -  Tl> 

JL0  ( I >  =  JL 
930  PRTT9U)  =  PPP 
K  s  0 

OC  940  I=1.MMASS 
II  =  ISP(I) 

OX  =  DSOC  <I) 

ppp  -  PPTTP(II)  *CX* (PRTTP  ( 1 1*1 )  -  PRTT0  (II)) 

940  PR  ( I )  =  PPP*ZZ 
C-C  TO  9000 

TAPE  OPTION. 

1000  IF  (NU.EC.O)  GO  TC  1 S 0 0 
C  UNIFORM  LOAD. 

CO  1100  K=JL. NTIME 
T 1  =  OEL^^LOAT (K-l) 

IF  (TIME.LE.T1)  C-C  TC  1200 
1100  CONTINUE 

JL  =  NTIME 

IF  (NTIME. LE.MAX02)  ppp  a  D ( l » NT IM£ ) 

IF  (NTIME. GT.MAXC2)  PPP  =  0 (1 ,NTIME-MAXC2) 

GO  TO  1300 


T1  a  (TINE-TD/CPL 
IP  ( Jl-1 .GT.MAXD2)  GC  TO  1220 
IF  ( JL.GT *maX02)  GO  TC  123o 
P2  =  P(l.JL) 

GC  TO  1250 

1220  ci  *  0 ( 1 ♦ JL-1-MAXC2) 

1230  P2  a  Q  ( l  ♦  JL-MAXC2) 

1250  PPP  a  pi  ♦  T 1  *  < P2-P 1 > 

1300  PPP  =  PPP*Z? 

IF  (KTYPE.LT. 6)  GC  TC  9000 

PX  =  PPP 

IF  (KTYPE.EO. 10)  FX  a  0. 

CO  1400  I*1.NVASS 
1400  PP (  I  )  a  PX 
GO  TO  9000 

C  NCN-UMFCfiM  LOAC. 

1500  CC  1600  KrJL.NTIYF 
T 1  a  OEL*PLOAT (K-l ) 

IF  (TIME.LE.T1)  GC  TC  1700 
1600  CONTINUE 

JL  a  NTIME 
GO  TO  1800 
1700  JL  =  K 

T 1  a  T 1  -OEL 
T 1  a  <TI*E-T1)/CEL 
1800  IF  (KTYPE.GT.5)  GC  TC  2400 
C  PANELS. 

IF  (JL.LT.NTIME)  GO  TC  I90o 
CO  1850  KGa  l.NC-SLM 

IF  (NTImE.LE.mAXC2>  PPTTB  (  kG )  =  P ( KG . NT  I ME ) 

IF  (NTIME.GT.MACX2)  PPTTB(wG)  r  0 ( KG »MT IMF-M A XD2 ) 
1850  CONTINUE 

GO  TO  2000 

1900  DO  1950  KGai.NGSL* 

IF  ( JL-1 .GT.maxC2)  GC  TO  1920 
PI  a  P (KG  * JL-1 ) 

IF  ( JL.GT .MAX02)  GC  TC  I93n 
P2  a  P(KC-tJL) 

GC  TO  1950 

1920  PI  a  0 (KG  » JL-1-m a  *02 ) 

1930  P2  =  0(KG»JL-MAXC2) 

1950  PPTTB  (KG)  a  p\  ♦  T1MP2-P1) 

C 

2000  K  s  0 

no  2300  1=1 ♦ ngt 
IF  (NPX.EG.l)  GC  TO  2050 
II  =  I  X  I  (  I  ) 

CX  a  0X1(1) 

2050  CO  2300  J=1,NPT 

IF  (NUSE(J.I) .EC.O)  CC  TC  ?300 


K  =  K  ♦  1 

IF  (NPY.EO.l)  GC  TO  2100 
JJ  =  JYJ(J) 

CYs  OY1 (J> 

PI  s  PRTTG ( JJ) 

ppp  =  PI  ♦  DY* (PPTTP ( JJ* 1 )  -  PI) 

PA(K)  =  PPP*7Z 
GC  TO  2300 

2100  IF  (J.GT.l)  GO  TO  2200 
P1=PRTTQ(II) 

PPP  s  (PI  ♦  OX<* (PRTT6  < I T ♦  1 )  -  Pl))*ZZ 
2200  P A ( K )  =  PPP 
2300  CONTINUE 

GO  TO  9000 
C  9EAMS • 

2^00  IF  < JL.LT.NTI*E)  GO  TC  2500 
no  2450  KG=1,NGSUM 

IF  (NTIME.GT.MAXC2)  FPtTP(k'G)  s  Q (KG,NTIME-^AXD2) 
2450  CONTINUE 

GO  TO  2600 

2500  DC  2550  KG=1,NGSU* 

IF  (JL-1.GT.MAXC2)  GC  TO  2c20 
PI  =  P(KG*JL-1) 

IF  (JL.GT.MAX02)  GO  TC  253n 
P2  =  P(KG.JL) 

GO  TO  2550 

2520  PI  =  Q<KGtJ|_-l-PAXC2> 

2530  P2  a  Q(KGtJL-MAX02) 

2550  PRTTB(KG)  a  PI  ♦  Tl*<P2-Pli 
2600  00  2700  I*1,NMASS 
II  =  ISP(I) 

OX  =  OSOC(I) 

PI  a  PRTTB(H) 

2700  PP(I)  a  ZZ*(P1  ♦  CX«<PRTTB(II*1)  -  Pll) 

C  * 

9000  RETURN 
END 
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APPENDIX  B 

RIB  BUCKLING  OPTION  CHANGES 

This  appendix  contains  the  changes  to  the  N0VA-2LT  (N0VA-2LTS)  codes 
which  are  necessary  to  include  pinned/sl iding-pinned  end  constraints  for  the 
rib  buckling  option  (KTYPE=10). 

Since  the  NOVA-2  and  N0VA-2S  program  libraries  are  maintained  at  the  AFWL 
using  the  Control  Data  Corporation  6000/7000  UPDATE  program,  the  UPDATE 
directives  needed  to  accomplish  these  changes  are  also  included. 
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1  i  iFr  Ifc 


"  nar  mi  iiaMMiB  .dap-  - 


CHaNC-FS  TC  EXTEKC 

N0VA2-LT  FOR  EXFER I^ENT A|  BUCKLING 


•IDEM  AUG0478 

♦INSERT  C0MP2.5I 

IF  (KEYB2.EQ.25  GC  TC  3970 
TH  =  ATAN2(W(l)-Wl.V(l)-Vl) 

STH2  =  SIN (TH) 

CTH2  =  CCS(TH) 

♦INSERT  COMP2.103 

IF  (KEYBl*KTYPE.EC.30>  IL1(I)  =  2 
♦INSERT  COMP2.114 

IF  (KEYB1.E0.3)  EL  =  P I/SCpT ( ( VM ( 1 ) -V2) »*2  ♦  (W* C 1 ) -W2) **2 ) 
♦INSERT  C0MP2.116 

IF  (KEY82.EQ.3)  IL  *  N 
♦INSERT  C0MP2.119 

IF  (KEYBl.EQ.2)  GC  TC  4075 

C  =  SORT ( ( VM (I )  -  VK<1))**2  ♦  ( WM  ( I >  -  WM  Cl)) **2 » 

0  =  AMP«SlN<D*EL> 

4075  CONTINUE 
♦INSERT  C0MP2.127 

SN(1)  =  .00298/CELTS <1) 

♦INSERT  CYCLE. 7 

IF  (KTYPE.EO.10)  PPP  =  PPP»WT*WR  ( 1 ) 

♦CELETE  CYCLE. 9 
♦INSERT  CYCLE. 52 

IF  (KTYPE.EQ.10)  CR I T  <  1 )  =  -1.0 
♦OELETE  EQUILP.S4 
250  CONTINUE 
♦INSERT  EQUILP.61 

IF  <KEY81*KTYPE,NE.3C >  RETURN 

C  PINNED-SLIDING  constraint. 

ACCNV(1)=C6(1)*(PFP*CCST(1 j ♦0IGN(2)*COST(2)-C(2)*SINT(2) ) 
ACCNW(l)  =  0.0 
•INSERT  FB .37 

KTYPEX  =  KJYPE 
KTYPE  =  7 
♦INSERT  FB .230 

KTYPE  =  KTYPEX 
•INSFRT  FB.236 

KTYPE  =  KTYPEX 
♦INSERT  FB . 242 

KTYPE  =  KTYPEX 
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CHANGES  7C  EXTE*C 

N0VA-2LTS  EOR  EXPERIMENTAL  BUCKLING 


•IOENT  AUG 0378 

•INSERT  C0MP2.S1 

IE  (KEYB2.E0.2)  GC  TC  3970 
TH  a  ATAN2(W{l)-Wl,V<l)-Vl) 

STH2  *  SIN ( TH> 

CTH2  =  CCS ( TH) 

•INSERT  COMP2.103 

IE  <KEY8l*KTYPE.EC.30)  IL  1  ( I )  a  2 
•INSERT  C0MP2.114 

IF  (KEYR1.E0.3)  EL  =  p I /SORT ( ( VM ( 1 ) -V2) **2  ♦ (WM ( 1 ) -W2) **2) 
•INSERT  C0MP2.116 

IE  (KEYB2.EQ.3)  IL  a  N 
•INSERT  COMP2.120 

IE  <KEYBl.EQ.2)  GC  TC  4075 

0  *  SORT ( ( VM ( I )  -  VM(1))**2  ♦  ( WM ( I )  -  WM ( 1 ) ) **2 ) 

0  a  AMP*SIN<0*EL) 

4075  CONTINUE 
•INSERT  COMP2.138 

SM ( 1 )  a  • 00339/CELTS  ( 1 ) 

•INSERT  CYCLE. 7 

IE  (KTYPE. EQ. 10)  PPP  *  Ppp*WT*WR (1) 

•DELETE  CYCLE. 9 
•INSERT  CYCLE. 52 

IE  (KTYPE.EQ. 10)  CRlT(l)  a  -l.o 
•INSERT  EQUILP.G1 

IE  (KEYB1*KTYPE.NE.3C>  RETURN 
C  PINNED-SLIDING  CONSTRAINT. 

ACCNV ( 1 ) aC6 ( 1 ) * (PPP«CCST ( 1 j ♦SIGN (2) *COST (2) -Q (2) *S INT (2) ) 
ACCNW(l)  a  o.O 
•INSERT  EB.37 

KTYPEX  a  KTYPE 
KTYPE  =  7 
•INSERT  EB.230 

KTYPE  =  KTYPEX 
•INSERT  EB.236 

KTYPE  a  KTYPEX 
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APPENDIX  C 


FREE  BOUNDARY  CONDITIONS  AND  DISCRETE  LINEAR  SPRINGS  OPTION  CHANGES 

This  appendix  contains  the  changes  to  the  NOVA-2LTS  code  which  are 
necessary  to  include  free  boundary  conditions  and  discrete  linear  springs. 

Since  the  N0VA-2S  program  library  is  maintained  at  the  AFWL  using  the 
Control  Data  Corporation  6000/7000  UPDATE  program,  the  UPDATE  directives  needed 
to  accomplish  these  changes  are  also  included. 


r>  n  r» 


EXTENO  ceprcp  to  include 
F-F  AND  C-F  EDGE  CCNDITIcNS  AND  ELASTIC  SPRINGS 


•IDENT  OCT0378 
•DELETE  JUNE2678.3 

3  PI  »X9  (28)  »XG  (28  )  ,N8N01  ,n6N02,NSPR »  IDIR  (30)  ,NSPC-  (30)  * 

4  NSPBOO)  ,8IGK(30)  «XtP3 
•DELETE  BOLT. <3 

DIMENSION  C01 (20) ,CC2  (20) ,cC3  (20) ,CD4(20) ,CD5(2C> ,CC6<20) 

•INSERT  BOLT. 20 

OATA  COS/O. 5P6864162ES8,  1.49417561428.  2.50024694616. 

1  3.49998931984.  4.50Q00046151 .  5.49999998001.  6.50000000087, 

6  7.49959999995.  8 .5 ,9.5 . 1 0 .5 . 1 1 .5, 12.5 » 1 3.5 ♦ 14 . 5 » 15.5 , 16 .5 , 

3  17. 5.18. 5. 19. S/ 

DATA  CO6/0. 734055513769.  1.01846731877,  0.999224496517, 

1  1.00003355325,  0.559998550107,  1.00000006264,  0.599999997254, 

2  1.00000000011,  12*1.0/ 

•INSERT  BOLT. 23 

FaCI  *  SCRT  (3.0) 

•DELETE  BOLT. 27 

GO  TO  (500,700.930.500,960),  NBND1 
•INSERT  BOLT. 65 

FREE-FREE*  GAMMA. 

930  OC  950  1*1, MG 
M  s  MGM(I) 

XI  =  COl(M) 

X2  *  CD2  (M ) 

CC  950  J=1,NGT 
II  *  II  ♦  1 
IF  (M-2)  935,940,545 
935  FP1  (II)  =  1.0 
FP2 (II)  *  0. 

FP3 ( 1 1 )  =  0. 

GC  TO  950 

940  FP 1  ( 1 1 )  *  FAC1*(1.0  -  2.0*pAM(J)/P1) 

FP2 (II)  *  -FAC  1*2. 0/E  1 
FP3(II)  =  0. 

GC  TO  950 

945  X3  *  X1*C-AM(J) 

F X 1  *  EXP ( X 3 ) 

EX2  a  EXP ( -X3 ) 

SL  *  SIN(X3) 

CL  *  COS(X3) 

FP 1  (II)  *  CL  -  X2*SL  ♦  .5«(1.-X2)*EX1  ♦  ,5»(1.  ♦  X2)*EX2 
FP2(II)  *  X 1  * ( -SL  -  X 2*CL  *  .5* ( 1 .-X2) *EX 1  -  ,5« ( 1 . ♦ X2 ) *EX2 ) 

FP3 ( 1 1 >  a  X 1 **2* ( -CL  ♦  X2*cl  ♦  .5* ( 1 . -X2 ) «EX 1  ♦  .5* ( l .♦XE) «EX2) 
950  CONTINUE 

CK (5)  a  1 . /FAC 
GO  TO  1000 
C 
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C  CLAmpfd-FRee,  Gamma. 

c 

960  DC  980  1=1, MG 
m  =  MGM(I) 

COL ( I )  =  CDS (M) 

<380  CCA  ( I )  =  CD6(M) 

GO  TO  540 
•DELETE  BOLT ,76 

GO  TO  (1100,1300, 153C, 1500, 1580) »  NRND2 
•INSERT  BOLT. 121 
C 

C  FREE-FREE,  BETA. 

C 

1530  00  1570  1=1, MB 
N  a  NBN(I) 

XI  =  CD1(N) 

X2  =  CD2(M 
DO  1570  J=1,NR1 
IT  =  II  ♦  1 

IF  (N-2)  1540, 1550, lSGo 
1540  FP5  (II)  =  1.0 
FP6  ( 1 1 )  =  0. 

F P 7  (II)  =  0. 

GO  TO  1570 

1550  FP5  (II)  =  F AC  1  * ( 1 .  -2 .*RETO ( J) /PI ) 

FP6 (II)  =  -FAC  1 *2 • 0/P 1 
FP7  (II)  =  0. 

GO  TO  1570 
1560  X2  a  X 1 *5ETR  ( J ) 

EX  1  a  EXP  <  X 3 ) 

EX2  a  EXP(-X3) 

SL  =  SIN ( X  3 ) 

CL  =  COS ( x  3) 

FP5 (II)  =  CL  -  X2«SL  ♦  ,5« ( 1.-X2) »EX1  ♦  ,5*(1 
FP6 (II)  =  X 1  * ( -SL  -  X2*CL  *  . 5» ( 1 . -X2) *EX 1  - 
FP7 ( 1 1 )  =  X1**2*(-CL  ♦  X2*-L  ♦  .5* ( 1 . -X2 ) *EX 1 
1570  CONTINUE 

CK(6)  =  1.0/FAC 
GO  TO  1600 
C 

C  CLAMPEO-FREE,  BETA. 

C 

1580  DC  1590  1*1, MB 
N  a  NBN  < I ) 

COL  ( I )  =  CD5  <N) 

1590  CO A ( I )  a  C06(N) 

GO  TO  1140 
•OELETE  DERV2.245 
C 

C  ELASTIC  SPRINGS. 

C 

IF  (NSPR.EQ.O)  GO  TC  1740 
00  1730  L 12  =  1,NSPR 
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♦  X2)»EX2 
5M1»*X2)*EX2> 

♦  , 5* ( 1 . ♦ X2 ) *EX2 ) 


1 


LI  a  N$PG  (L12) 

L2  =  NSPf?(L12) 

K  =  (L1-1)*NRT  ♦  L 2 
8TGKL  *  eiGK<Ll2) *XLP3 
IF  ( IDIR (L 12)  -  2)  1710,1720.1725 
1710  SUPS  =  BlGKL*U(K)«SlNGtHHG.Ll)*SIK2B(NK,0*L2)  ♦  SUPS 
GO  TO  1730 

1720  SVRS  =  8IGKL»V(K)*SlN2G(MMo/*Ll)*SINB(NNO*L2)  ♦  SVRS 
GO  TO  1730 

1725  SWPS  =  BIGKL*W(K)«FPl <MMC*Ll)»FP5(NNO*L2)  ♦  SWR$ 

1730  CONTINUE 

1740  IF  (ARS(SWRS) .GT.1.E20)  GO  TO  2150 
•DELETE  DSET1.53 

190  READ  (5.7000)  NSPR 

IF  (NSPR.EQ.O)  GO  TC  210 
DO  200  1  =  1 , NSPP 

200  READ  (5.7110)  IC  IR  ( I  > *NSPG  ( I ) .NSPR ( I ) ,B IGK ( I ) 

210  KSTIF  a  IABS(NSG)  ♦  IABS(NSS) 

•DELETE  0SET1.164 
1190  WRITE  (6,12000)  NSPR 

IF  (NSPR.EQ.O)  GO  TC  1210 
WRITE  (6.12200) 

00  1200  1=1, NSPR 

1200  WRITE  (6.12100)  IC IP <1 > ,NSPG ( I > .NSP8 < I ) ,B IGK ( I) 

1210  IF  (NSG.EQ.O)  GO  TC  1300 
•INSERT  OSET1.219 

NRNOl  a  NBNO/IO 
•INSERT  0SET1.223 
7110  FORMAT  (3I12.F12.1) 

•DELETE  0SET1.225 

7170  FORMAT  (75H1 INPUT  OATA  FOR  CERROP  (MODIFIED  TO  INCLUDE  EXTRA  B.C 
1  1  AND  ELASTIC  SPRINGS) > 

i  *  INSERT  0SET1.276 

12000  FORMAT  (8H0NSPR  *  13) 

[  12100  FORMAT  (5X, 316, E15.6) 

12200  FORMAT  (37h  IOlR 

•OELETE  DSET3.32.DSF13.37 
>  2950  IF  (NRNOl. EQ.l)  WRITE 

IF  (N8N01.EQ.2)  WRITE 
IF  (NBN01.EQ.3)  WRITE 
IF  (NRNOl. EO. 4)  WRITE 
i  |  IF  (NBN01.E0.5)  WRITE 

*  IF  (NBN02.E0.1)  WRITE 

IF  (NRN02.E0.2)  WRITE 
IF  (N8N02.E0.3)  WRITE 
IF  (NRN02.E0.4)  WRITE 
IF  (NBN02.E0.5)  WRITE 
•INSERT  0ESET3.73 

WRITE  (6*13400)  NSPR 
•OELETE  DSET3.ee, 0SET3.E9 

IE  (NRNOl. GT. 3)  NNSY^G  a  1 
IF  (NBND2.GT.3)  NNSYMB  a  1 
•INSERT  0SET3.126 
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(6,9r90) 

(6.10000) 
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IF  (NSPR.EO.O)  GC  TC  ^075 
XLP3  =  2.0»XJ/<A*FI*H*XLP> 

IF  <NSYMG*NSYNB.EC.1>  GO  TO  4075 


ACCOUNT  FOR  SYMMETRY* 

XLP3  *  4 •  0*XLP3/  FLCAT(NSY(-e*l)*(vSYNBM)) 

CC  4060  L=1«NSPR 
I  =  NSPG(L) 

j  s  NSPB  <L)  _ 

IF ( I .EQ.mBaP)  PICK  ( L )  =  0,c«BIGK (L) *FLCAT (NSYMG*1 > 

IF  (J.EQ.NSAR)  0  TGK  (L  >  =  0 .5«RIGK (L) *ELOAT <NSYMP*1 ) 

4060  CONTINUE 
4075  CONTINUE 

•CELETE  MAY0878 • 31  _  , 

9300  FORMAT  (1M.20X*33HC  E  P  R  C  P  (EXTENOEC  CCT.«  1978)/ 

1  15H0PANEL  ANALYZED 

•INSERT  DSET3.231  _ 

9930  FORMAT  (31h  FREE  •  FREE*  GAMMA  DIRECTION) 


9930  FORMAT  (31h  FREE  •  FREE*  GAMMA  DIRECTION) 

•  INSERT  OSET1.232  „ 

9950  FORMAT  (34h  CLAMPEC  -  FRFE  ♦  GAMMA  DIRECTION) 

999flPF0RMAT  1 30h  FREF  -  FREE*  BETA  DIRECTION) 

10050  FORMAT  (23H  CLAMPEC  -  FRFE*  BETA  DIRECTION) 

1 34 00^ FORMAT^ *( 3BH0NUMEER  CF  DISCRETE  ELASTIC  SPRINGS  *  13) 
•DELETE  DTSTEP.19*0TSTEP.2C 

IF  (NRN01.EG.4)  Cm  =  *15 
IF  (N0ND1 *E0« 1 )  CM  =  .3 
•OELETE  DTSTEP.22,DTSTER,23 

IF  (NBN02.E0.4)  CN  =  *15 
IF  (NRND2.E0.1)  CN  =  *3 
•IOENT  AUG1678 

5TP1  (10*2)  *STRA  (5*21*20)  »ST1*ST2*TEN01*TEND2*TETA  (11)  * 

*10  CCT0278 
♦I  PIMT.54 


CLAMPEC  -  FRFE*  GAMMA  DIRECTION) 
FREF  -  FREE*  BETA  DIRECTION) 
CLAMPEC  -  FRFE*  BETA  DIRECTION) 
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